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Abstract 

Recent developments in the evaluation of vector-meson spectral functions in hot 
and dense matter are discussed with emphasis on connections to the chiral phase tran- 
sition in QCD. Model independent approaches including chiral low-density expansions, 
lattice QCD, chiral and QCD sum rules are put into context with model predictions for 
in-medium vector-spectral function utilizing effective Lagrangians. Hadronic many- 
body calculations predict a strong broadening (and little mass shift) of the p spectral 
function which rapidly increases close to the expected phase boundary of hadronic 
and quark-gluon matter. Pertinent dilepton rates appear to degenerate with pertur- 
bative quark-antiquark annihilation in the Quark-Gluon Plasma, suggestive for chiral 
symmetry restoration. Applications to low-mass dilepton spectra in heavy-ion colli- 
sions result in quantitative agreement with recent high-quality data at the CERN-SPS. 
Thermal radiation from temperatures around T c consistently reproduces the experi- 
mental dilepton excess observed at masses above 1 GeV as well. The interpretation 
of dilepton sources at high transverse momentum appears to be more involved. 
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1 Introduction 



The quest for the elementary entities of matter has always been a central objective in 
physics. Of no less interest is the emergence of the structure (or phases) of matter built 
from its basic constituents and their interactions. The exploration of matter governed 
by the strong force is at the forefront of contemporary research in nuclear physics. The 
most common form of strongly interacting matter in the present-day universe is contained 
in atomic nuclei, which are bound states of nucleons, i.e., protons and neutrons. Nuclear 
matter as found in the center of heavy nuclei is characterized by a nucleon (energy) density 
of about qq = 0.16 fm~ 3 (eo = 0.15 GeV/fm 3 ), rendering one table spoon of this material 
a mass of about one million kilotons (10 12 kg). The binding of nucleons is strong enough 
to cause a reduction of nuclear masses by about 1% compared to the sum of the rest mass 
of the individual nucleons, — 940MeV/c 2 . But how does the mass of a nucleon arise? 
In the late 1960's it was discovered that the nucleon itself is a composite object, built of 
three "valence" quarks of up (u) and down (d) "flavor" . The bare masses of u and d quarks 
are only about 5-10MeV/c 2 , and believed to be generated by a condensate of (yet to be 
discovered) Higgs bosons in the electroweak (EW) sector of the Standard Model of Ele- 
mentary Particles. That is, about 98% of the nucleon's mass is generated dynamically by 
the strong interaction. Moreover, no individual quarks have been observed in nature thus 
far: they are "confined" into hadrons, either baryons or mesons (conglomerates of three 
valence quarks or of a quark and antiquark, respectively). In the 1970's, the quantum 
field theory underlying the strong force has been developed, Quantum Chromodynamics 
(QCD), based on quarks and gluons as fundamental degrees of freedom. This theory has 
been quantitatively confirmed in high-energy scattering experiments, where the strong 
coupling constant, a s ~ 0.1, is relatively small and perturbation theory can be reliably 
utilized to obtain quantitative results for observables. However, at low momentum trans- 
fers, a s becomes large, perturbation theory ceases to be applicable and nonperturbative 
mechanisms take over. It is in this regime where quark confinement and mass generation 
occur, posing formidable challenges for their theoretical understanding PQ. 

1.1 QCD Vacuum and Chiral Restoration 

It turns out that confinement and mass generation are intimately connected with the phase 
structure of strongly interacting matter (see, e.g., Ref. [2]). In fact, even the structure of 
the QCD vacuum is far from trivial: similar to the EW sector, it is believed to be filled 
with condensates, which are closely related to the origin of hadronic masses. There are, 
however, important differences: the QCD condensates are made of (scalar) composites of 
quarks and gluons (rather than elementary fields like the Higgs boson), and they do not 
induce a breaking of the gauge symmetry. In what follows, the so-called "chiral" quark- 
antiquark condensate, (Olf/gjO) ~ (— 250 MeV) 3 , will be of particular importance. It breaks 
the (approximate) chiral symmetry of QCD, which corresponds to the conservation of left 
and right "handedness" of massless quarks (applicable for the light u and d quarks, whose 
masses are parametrically small, m u ^ <C | (0|g(/|0) j 1 / 3 ). While the quark condensate cannot 
be directly observed, its consequences are apparent in the excitations of the condensed 
ground state, i.e., in the hadron spectrum. Since chiral symmetry is a global symmetry 
(rather than a local one depending on space-time position), its spontaneous breaking must 
be accompanied by (almost) massless Goldstone bosons. For two quark flavors the latter 
are identified with the three charge states of the pion, whose mass, m n ~ 140 MeV, is 
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"abnormally" small compared to that of all the other hadrons (e.g., ~ 940 MeV). 
The observed hadron spectrum encodes further evidences for the spontaneous breaking 
of chiral symmetry (SB^S): chiral multiplets (e.g., p(770)-ai(1260) or iV(940)-7V*(1535)), 
which would be degenerate if the ground state were chirally symmetric, exhibit a large 
mass splitting of typically AM ~ 500 MeV. The effects of SB^S seem to (gradually) cease 
as one goes up in mass in the hadronic spectrum [3]. This is one of the indications that 
SB^S is a low-energy, strong-coupling phenomenon which is no longer operative at high 
momentum transfers where perturbation theory becomes applicable. 

When heating the QCD vacuum its condensate structure is expected to change. Loosely 
speaking, thermally excited hadrons "evaporate" condensed qq pairs which eventually 
leads to the restoration of the spontaneously broken chiral symmetry. Numerical compu- 
tations of the lattice-discretized path integral for QCD at finite temperature predict chiral 
symmetry restoration (xSR) to occur at a (pseudo-) critical temperature of T c ~ 160- 
190 MeV [HE]; corresponding to an energy density of about e c ~ lGeV/fm 3 . The chiral 
transition is characterized by a rapid decrease of the qq condensate, which, in fact, serves 
as an order parameter of strongly interacting matter. In the limit of vanishing light quark 
masses and for three quark flavors, this transition is of first order, while for realistic quark 
masses as realized in nature (two light quarks u and d and a more heavy strange quark, 
m s ~ 120 MeV), it is more likely a rapid cross-over. Key manifestations of chiral symme- 
try restoration are its (observable) consequences for the hadron spectrum. Chiral partners 
must degenerate implying massive medium modifications of hadronic spectral functions 
as the transition is approached. This notion is a quite general concept found, e.g., in 
solid state physics where phase transitions are routinely diagnosed utilizing "soft-mode 
spectroscopy". This applies in particular to a second order phase transition where the 
mode associated with an order parameter becomes massless (soft). But even for bulk 
matter properties, rapid changes in the thermodynamic state variables are directly related 
to changes in the relevant degrees of freedom at the typical thermal scale (temperature or 
Fermi momentum). Interestingly, the chiral transition is accompanied by the dissolution 
of hadrons into quarks, i.e., the deconfinement transition, at the same temperature (at 
least for vanishing net baryon density) . The reason for the apparent coincidence of the two 
transitions is not understood. The deconfined and chirally restored strongly interacting 
matter is commonly referred to as the Quark-Gluon Plasma (QGP). The experimental ver- 
ification and theoretical understanding of the mechanisms leading to the QGP are central 
objectives in modern nuclear research. 

1.2 Ultrarelativistic Heavy-Ion Collisions 

The only way to produce and study hot and dense strongly interacting matter in the 
laboratory is by colliding atomic nuclei at high energies. Several large-scale experiments 
at ultrarelativistic bombarding energies, 3> m^, have been conducted over the past 
~20 years, most recently at the SPS at CERN (at center-of-mass energies up to \/s = 
17.3 AGeV) and at the Relativistic Heavy-Ion Collider (RHIC) at Brookhaven (up to 
y/s = 200 AGeV) [7] with the heaviest available nuclei at A~200 (Pb and Au). 

The first question that needs to be answered is whether these reactions produce equili- 
brated matter, i.e., do the produced particles undergo sufficient rescattering to justify the 
notion of an interacting medium characterized by bulk thermodynamic variables? Exten- 
sive and systematic measurements of hadronic observables have lead to a positive answer. 
This is extremely exciting as it puts within grasp the possibility to recreate, at least for a 
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short moment, the matter which the early universe was made of just a few microseconds 
after the Big Bang! While hadronic measurements are discussed and interpreted in depth 
in other contributions of this volume, let us sketch some of their main features: transverse- 
momentum (pt) spectra of different hadron species (pions, kaons, protons, etc.), which 
characterize the hadronic fireball just before break-up at its "thermal freeze-out", exhibit 
a collective explosion reaching an average speed of about half the speed of light at a final 
temperature of T f ^ h ~ 100 MeV. The ratios of the observed hadron species point at a sig- 
nificantly higher temperature of T^ 1 ~ 160 MeV [8j [9], i.e., the chemistry of the fireball 
(driven by inelastic scattering processes) appears to freeze out significantly earlier than 
kinetic equilibrium (maintained by elastic interactions). This is consistent with the large 
difference of empirical elastic (e.g., irir — > p — > irir or ttN — > A — ► ttN) and inelastic (e.g., 
7T7T — > KK) hadronic cross sections, with typical values of ~100mb vs. ~lmb, respec- 
tively. Since the cross sections determine the relaxation times according to r ~ {Qh^Vrd) 1 
(gh- hadron density, v rc y. relative velocity of the colliding hadrons), one obtains a clear 
hierarchy in the underlying relaxation times, r t h <C r^. The interacting hadronic phase 
between chemical and thermal freeze-out will play an important role in the remainder of 
this article. More differential analyses of the flow patterns of the measured hadrons allow 
to trace back the matter properties to earlier times in the evolution of the fireball. In par- 
ticular, the magnitude of the "elliptic flow" measured at RHIC indicates that the medium 
thermalizes on a rather short time scale, r < 1-2 fm/c after initial impact, translating into 
(energy-) densities of a factor 10 or more above the critical on^. A thermal (hydrody- 
namic) description [1Q|. \TT\ [I2j [13] of the fireball in semi- / central collisions of heavy nuclei 
at RHIC appears to be valid for hadrons up to momenta of pt — 2-3 GeV, comprising 
approximately 95% of all produced particles. At high transverse momenta, px > 5 GeV, 
hadron production is dominated by hard scattering, i.e., a primordial parton-parton colli- 
sion at high momentum transfer within the incoming nucleons, followed by fragmentation 
into (a spray of) hadrons (jets). In central Au-Au collisions at RHIC, a factor of ~5 
suppression of high-pr hadron production has been observed ("jet quenching" j^. While 
these hadrons (or their parent quarks) do not thermalize, their suppression indicates a sub- 
stantial coupling to the created medium, associated with an energy loss of a fast parton 
propagating through the fireball. The (energy-) density of the medium required to ac- 
count for this effect is roughly consistent with the estimate inferred from a hydrodynamic 
description of the elliptic flow of Iow-pt hadrons. 

A second level of questions concerns the relevant degrees of freedom of the produced 
matter, i.e., whether there is explicit evidence that individual partons leave a distinctive 
footprint in observables. It turns out that the elliptic flow is again revealing interesting 
features in this context: it has been found |144 I15| that the elliptic-flow coefficient, v^Kt), 
of all measured hadrons, h=ir, K, p, A, S, <fi, exhibits a remarkable universality as a 
function of transverse kinetic energy, Kt=vtit — (mT=(pj^+mf l ) 1 ^ 2 ): when scaled with 

1 Elliptic flow characterizes the azimuthal asymmetry in the p-r-spectra of particles (in the plane trans- 
verse to the beam axis). In a non-central heavy-ion collision, the initial nuclear overlap (interaction) zone 
is "almond-shaped". If the system thermalizes before this spatial anisotropy is smeared out (e.g., due to 
free streaming), a larger pressure gradient builds up along the "short"' compared to the "long" axis of 
the initial almond. This thermal pressure drives a collective expansion of the "almond" which is stronger 
along the short axis and thus results in particle momenta with a preference to be aligned with this axis. 
The magnitude of the elliptic flow is thus sensitive to how fast thermalization is established. 

2 Jet quenching is probably also present at the SPS but it is quantitatively smaller than at RHIC (about 
a factor of 2 suppression) and masked by a large initial pr broadening in the interpenetrating nuclei prior 
to the hard scattering, known as "Cronin effect". 
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the constituent-quark number content, n q , of hadron h, all measured hadron-f2 data ap- 
pear to collapse on a single curve, v^K^Kt / n q )=V2 (Kt) /n g . This has been interpreted 
as evidence for a collectively expanding partonic source hadronizing via quark coalescence. 
A fully consistent theoretical description of this phenomenon has not been achieved yet. 

A third level of investigations has to address signals of the deconfinement and/or 
chiral restoration transitions. In a rigorous sense, this requires the assessment of order 
parameters associated with these transitions. However, changes in order parameters are 
not always easily observable. This is particularly true in the present context and we are 
led back to the idea of "mode spectroscopy", to be conducted in the environment of a 
short-lived, rapidly expanding fireball of a heavy-ion collision. Individual (stable) hadrons 
emanating from the collision zone have all long recovered their free (vacuum) masses by 
the time they are measured in the detectors. A better observable are invariant-mass 
spectra of short-lived resonance decays, h — ► h±h2, with a lifetime, t/j, comparable to, or 
smaller than, the lifetime of the interacting fireball, tfb — 10 fm/c. Such a resonance (e.g., 
A — ► ttN or p — > irir) has a large probability to decay inside the medium so that its decay 
products can carry the information on its invariant- mass, m? = (pi +P2) 2 , at the point 
of decay to the detector. In principle, this would allow to determine the invariant-mass 
distribution (or spectral function) of the resonance h in the medium. The problem is that 
the decays products, h\ and h,2, are likely to undergo further rescattering in the fireball 
which destroys the desired invariant-mass information. The latter will thus be largely 
restricted to the dilute (break-up) stages of the medium in a heavy-ion collision. 

1.3 Dilepton Spectroscopy 

The decisive step to obtain access to hadronic spectral functions in the hot and dense 
regions of the medium is provided by electromagnetic (EM) probes, i.e., photons (7) and 
dileptons (arising from virtual (timelike) photons, 7* — ► l + l~ with l=e or ji) |16 (. 117 ^ [T8]. 
These are not subject to the strong force and thus suffer negligible final-state interac- 
tions, with a mean free path which is much larger than the typical size of the fireball, 
i?FB — 10 fm. The natural candidates for in-medium spectroscopy are the vector mesons 
(V), which carry the quantum numbers of the photon (spin-parity J p =l~) and thus 
directly couple to exclusive dilepton final states, V — > l + l~. In the low-mass region 
(M<1 GeV), which is the region of interest to study chiral restoration, the prominent vec- 
tor mesons are p(770), cj(782) and (/>(1020). In fact, the famous vector dominance model 
(VDM) pjj] asserts that the coupling of a (real or virtual) photon to any EM hadronic 
current exclusively proceeds via an intermediate vector meson (which is excellently sat- 
isfied in the mesonic sector but subject to corrections in the baryonic sector). Thus, 
if VDM holds in hadronic matter, dilepton emission is indeed equivalent to in-medium 
vector-meson spectroscopy. In thermal equilibrium, the contribution from the (isovector) 
p meson dominates over the u by a factor of ~10 (factor ~5 over the <f>, which, however, is 
thermally suppressed due to its larger mass). Furthermore, it can be shown [20J that, in 
the context of a heavy-ion collision, low-mass dilepton radiation from the hadronic phase 
dominates over the emission from a putative QGP phase, even at collider energies (RHIC 
and LHCjl. The excitement (and theoretical activity) in the field was further spurred by 
the suggestion of Brown and Rho [21] that the p-meson mass should drop to (almost) zero 

3 This is due to the much larger three- volume in the hadronic phase; at larger dilepton masses, M>1 GeV, 
the thermal Boltzmann factor, e~ 90/ ' T ', augments the sensitivity to higher temperatures which increases 
the QGP contribution relative to the hadronic one in the inclusive dilepton spectrum. 
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Figure 1: Left panel: ratio of cross sections, R = CTee^hadrons/oee-*/^; for electron-positron 
annihilation into hadrons relative to muon-antimuons, as a function of center- of-mass 
energy, y/s. The experimental data exhibit a nonperturbative resonance regime up to 
\/s~ 1.1 GeV, exhausted by the light vector mesons p, oj and (f), followed by a transition to 
an almost structureless perturbative regime at \fs>1.5 GeV. The latter is well described by 
perturbative QCD (pQCD), especially for \fs>2 GeV, where residual "oscillations" (due 
to excited vector resonances) have essentially ceased. The naive quark-model prediction 
(leading- order pQCD, 0{cP s )) is given by Rqm =N c ^2 q=u d s e 2 q =2. Right panel: strong 
coupling constant as a function of momentum transfer I33\j : note the increase of a s toward 
small Q suggestive for the emergence of nonperturbative phenomena. 



as a consequence of xSR- Early dilepton measurements in S(200 AGeV)-Au collisions at the 
CERN-SPS by the CERES collaboration [22] found a large enhancement of the spectrum 
at invariant masses below the nominal p mass, i.e., for M~0.2-0.7 GeV. These data could 
be well described by a dropping-mass scenario implemented into relativistic transport mod- 
els within a mean-field description [23\ [24"] . Subsequently, more "conventional" medium 
modifications of the p meson were investigated based on its rescattering on constituents 
of a hadronic medium, see, e.g., Refs. |25[ 126 1 I28j for reviews. The generic finding of 
these hadronic many-body calculations was a strong broadening of the p spectral function, 
which, when extrapolated to the putative phase transition temperature, T c ~ 175 MeV, 
leads to a complete "melting" of the resonance structure [29J . The broadening effect could 
account for a large part of the low- mass dilepton excess observed in S-Au collisions. The 
agreement was even better |30j with improved CERES/NA45 measurements carried out 
in the heavier Pb(158 AGeV)-Au system [31] 132] . The connection of the p melting to xSR 
appeared to be less direct than in dropping-mass scenarios. However, in Refs. |34[ [30] 
it was found that the hadronic dilepton rates following from the "melted" p close to T c 
rather closely resemble the rates computed in a partonic description, i.e., perturbative 
quark-antiquark annihilation. In the vacuum, such a phenomenon is well known from the 
e + e~ annihilation into hadrons: for M>1.5GeV the total cross section is well described 
within perturbative QCD using quark-antiquark final states, known as "parton-hadron 
duality" , cf. Fig. [1] It was therefore suggested that the conceptual implication of the p 
melting is a reduction of the "duality threshold" [Ml GO], from Md ua i=l-5 GeV in the vac- 
uum to essentially zero around T c . Note that a "perturbative" dilepton rate automatically 
implies chiral restoration (i.e., degeneracy of vector and axialvector channels). 

The accuracy in the 1995/1996 CERES /NA45 dielectron data [Ml E2] did not allow 
for a decisive experimental discrimination of the dropping-mass and melting-resonance 
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scenarios. An important step forward was realized with the NA60 dimuon spectra [35] in 
In(158 AGeV)-In collisions at the SPS. Excellent mass resolution and superior statistics 
enabled, for the first time, an isolation of the "excess radiation" (by subtraction of final- 
state hadron decays). The shape of the excess spectrum clearly favors a broadened p 
spectral function over scenarios involving dropping masses. The original predictions of 
hadronic many-body theory [30|, [20] are, in fact, in quantitative agreement [36] with the 
inclusive mass spectra in semi-/central In-In collisions. In the last round of CERES/NA45 
data [37] excess spectra have also been extracted in Pb-Au collisions (by subtraction of 
final-state hadron decays using a statistical model (Sj). While the overall data quality 
does not reach the level of NA60, the larger collision system and the access to very small 
dilepton masses in the dielectron channel (dimuons have a threshold of 2771^=210 MeV) 
can provide additional insights. 

The dilepton program at the CERN-SPS has thus far reached the highest level of 
maturity in the heavy-ion context. It also included a CERES/NA45 measurement in 
a low-energy Pb(40 AGeV)-Au run [38], which produced tantalizing hints for an even 
larger excess than at 158AGeV, but was unfortunately hampered by low statistics. At 
much lower, relativistic bombarding energies (l-2AGeV), the DLS collaboration at the 
BEVALAC reported a very large dilepton excess [39], which has recently been confirmed by 
the HADES collaboration at SIS [A0\ I41j. On the other hand, the dilepton measurements 
at RHIC are still in their infancy (first data indicate substantial excess radiation [42]), 
but it will become a central component in future runs [33]. Very interesting results are 
also emerging from vector-meson spectroscopy in cold nuclei using elementary projectiles, 
i.e., photons [44 \ 143] or protons [16]. It turns out that all of these observables are closely 
related, and their broad understanding is essential for the determination of the in- medium 
vector-meson spectral functions. Of particular importance is the consistency of theoretical 
descriptions beyond phenomenological applications and the interrelations between different 
approaches (including effective hadronic and quark models, lattice QCD and constraints 
from sum rules), which will ultimately reveal the mechanisms of chiral restoration. In this 
article, we give an up-to-date account of these efforts with special emphasis on a broader 
picture in the context of xSR- 

1.4 Outline 

Our article is organized as follows. In Sec. El we start by recollecting basic features of 
spontaneous chiral symmetry breaking in the QCD vacuum with emphasis on condensate 
structures and consequences for the hadronic excitation spectrum (sub-Sec. 12. ip . followed 
by a discussion of in-medium condensates within the landscape of the QCD phase dia- 
gram (sub-Sec. I2.2p . In Sec. El we scrutinize the links of the chain with which one hopes 
to connect thermal dilepton rates and (partial) xSR- We first introduce the EM correlation 
function which is the basic quantity figuring into the thermal dilepton rate (sub-Sec. I3.ip . 
Model-independent evaluations of medium effects can be obtained in the low-density limit 
from current algebra, in the high-temperature limit from perturbative QCD and, for van- 
ishing baryon-chemical potential from lattice QCD (sub-Sec. I3.2() . A valuable source of 
model-independent constraints is provided by chiral and QCD sum rules (sub-Sec. 13, 3[) 
which are energy moments of spectral functions that directly relate to order parameters of 
QCD and are generally not restricted in temperature and density. For practical applica- 
tions, effective hadronic models are an indispensable tool (sub-Sec. I3.4j) : their reliability, 
based on the choice of interaction vertices and associated parameters, crucially hinges 
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on a thorough procedure of theoretical and phenomenological constraints; an important 
question will also be the fate of the vector dominance model in the medium. In Sec. 0] the 
theoretical developments are tested in recent dilepton production experiments, starting 
with elementary reactions off nuclei representative for medium effects in cold nuclear mat- 
ter (sub-Sec. 14. The main part of Sec. U]is devoted to an analysis of dilepton spectra 
in ultrarelativistic heavy-ion reactions (sub-Sec. 14. 2p . focusing on recent results obtained 
at the CERN-SPS by the NA60 and CERES /NA45 collaborations. The spectral analysis 
is completed by a critical assessment of the combined theoretical and experimental status 
to date (sub-Sec. I4.3j) . We finish with concluding remarks in Sec. 

2 Chiral Symmetry, Condensates and Chiral Restoration 

It is generally accepted that strong interactions are described by Quantum Chromody- 
namics (QCD), introduced in 1973 [471 SHI ES] , with a Lagrangian density given by 

Cqcd = q(i-fD^ - M q )q - \g%G^ , D, = d, + ig s ^A« , (1) 

formulated in terms of elementary quark (q) and gluon {A^) fields (7^ and A a : Dirac and 
Gell-Mann matrices, respectively, A^„=diag(m u , m^, . . . ): current-quark mass matrix). In 
addition to the local SU(3) color gauge symmetry, Cqcd possesses several global symme- 
tries. The most relevant one in the present context is Chiral Symmetry, which can be ex- 
hibited by rewriting Cqcd in terms of left- and right-handed quark fields, qL,R=^{^- z fj5)q' 

£ Q CD = qLi^D^L + q-Ri^D^R - (q L M q q R + q R M q q L ) - ^G^G^ . (2) 

For small quark masses, i.e., u and d quarks, £qcd is approximately invariant under 
rotations qi,R — > e~ iaL ' R ' T ^ 2 qL,R, where or,l are 3 real angles and r operates in (u-d) 
isospin space. Chiral invariance of the QCD Lagrangian thus refers to the conservation 
of quark handed-ness and isospin. Alternatively, one can rewrite the chiral rotations 
as q — ► e~ lo ' v " r / 2 q and q — > e~ l ' y5aA " r ^ 2 q, giving rise to conserved isovector-vector and 
-axialvector currents, 

Jv = on^i . T A = qi^i^q ■ (3) 

2.1 Condensates and Hadron Spectrum in Vacuum 

As emphasized in the Introduction, the nonperturbative structure of the QCD vacuum 
is characterized by its condensates. A special role is played by the quark-antiquark (qq) 
and gluon (G 2 ) condensates. Apart from being the condensates involving the minimal 
number of quark- and gluon-fields, the former is a main order parameter of SB^S while 
the latter dominantly figures into the energy-momentum tensor of the theory. The vacuum 
expectation value of the latter's trace is given by 

(Tp=e-3P=-(G 2 ) + m q (qq) (4) 

where G 2 = —(/3(g s )/2g s ) Ga V G a ^ u involves the gluon-field strength tensor and the renorma- 
lization-group beta function, f3(g s ). The latter appears because the nonvanishing vacuum 
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Figure 2: Experimental data \50$ of the isovector-vector (IJ P =1\~, left panel) and 
isovector-axialvector (I J p =11 + , right panel) spectral functions from hadronic decays of 
t leptons (produced in y/s=91 GeV e + e~ annihilation at LEP) into even and odd numbers 
of pions/kaons, respectively. The lines indicate theoretical calculations using pQCD. 

value of Tjf breaks the scale invariance of the classical QCD Lagrangian, induced by quan- 
tum loop corrections. The small current light-quark masses, m q ~5 MeV, render the con- 
tribution of the quark condensate to Tjf small. The absolute value of the gluon condensate 
is not precisely known, but presumably rather large, around 1.5 GeV/fm 3 ~(330 MeV) . 
In fact, the magnitude of the quark condensate is not small either, (gg)~(-250 MeV) 3 
per light-quark flavor, and about 50% of that for strange quarks. This implies that the 
vacuum is filled with ~5 quark-antiquark pairs per fm 3 ! Also note that the quark con- 
densate maximally violates chiral symmetry by mixing right- and left-handed quarks, 
{qq)={qlQr + QrQl), implying that a quark propagating through the vacuum can flip its 
chirality by coupling to the condensate. The intimate relation between chiral symmetry 
breaking and the associated Goldstone-boson nature of the pion is highlighted by the 
Gell-Mann-Oakes Renner (GOR) relation, 

m lfl = -2m q (qq) , (5) 

which combines the effects of explicit chiral symmetry breaking, m^ocrnq, and SB^S with 
the pion decay constant as order parameter. 

One of the best direct empirical evidences for the spontaneous breaking of chiral 
symmetry is found in the vector channel, more specifically the isovector-vector channel 
(IJ P =11 + ) and its chiral partner, the isovector-axialvector one (IJ P =11~), precisely the 
Noether currents in Eqs. ([3j) . The pertinent spectral functions have been measured with 
excellent precision (and a detailed decomposition of the hadronic final states) at the Large 
Electron-Positron collider (LEP) in hadronic r decays by the ALEPH [50] and OPAL [51] 
collaborations, cf. Fig. [2l In the low-mass region, the strength of each of the two spectral 
functions is largely concentrated in a prominent resonance, i.e., the p(770) and ai(1260). 
This very fact indicates that the low-mass regime is dominated by nonperturbative ef- 
fects, while the (large) difference in mass and width of these resonances signals chiral 
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Figure 3: Schematic QCD phase diagram including empirical extractions of (// b,T) -values 
from observed hadron production ratios in heavy-ion experiments at different beam ener- 
gies J2|/; the bands indicate lattice- QCD and model estimates of the transition regions be- 
tween HM, QGP and CSC; with the HM-QGP transition (along the finite T, hb=0 axis) 
presumably a cross-over, and the HM-CSC one (along the finite [is, T=0 axis) possibly first 
order, there is presumably a second order endpoint, e.g., around (/fB,T)~ (400,160) MeV. 



symmetry breaking. This connection can be quantified by chiral sum rules developed by 
Weinberg [52] and others [53] in the late 1960's based on current algebra of chiral sym- 
metry. These sum rules relate moments of the difference between vector and axialvector 
spectral functions to chiral order parameters. In the chiral limit (777.^=0) one has 

oo 

r ds 

fn = ~ —s n [Imn y (s) - Imll^O)] , (6) 



f_ 2 = fl^L-F A , f-x = fl, / = 0, f 1 = -2ira s (0) (7) 

(r w : pion charge radius, F A : coupling constant for the radiative pion decay ir^ — > /j ± i/ m 7, 
four-quark condensate). 



2.2 Phase Diagram and Chiral Restoration 

A schematic view of the QCD phase diagram is displayed in Fig. [3l It is roughly character- 
ized by three major regimes (all of which most likely exhibit rich substructures): hadronic 
matter (HM) at small and moderate temperature (T) and baryon chemical potential (/is), 
Quark-Gluon Plasma (QGP) at high T and Color Super-Conductors (CSCs) at high /jg 
but low T. The latter may occur in the core of neutron stars, but are unlikely to be 
produced in heavy-ion collisions and will not be further discussed here. 

A key issue toward understanding the phase structure of QCD matter is the tem- 
perature and density dependence of its condensates. Various condensates serve as order 
parameters of broken symmetries and govern the (hadronic) excitation spectrum. The lat- 
ter provides the connections to observables. A first estimate of the medium modifications 
of the condensates can be obtained in the low-density limit [551 [56] , by approximating the 
thermal medium by non- interacting light hadrons, i.e., pions at finite T an nucleons at 
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finite qn. For the quark condensate, this leads to a linear density expansion of the type 



((qq))(T,HB) V QjXh „, T 2 1 Q N 

where denotes the scalar density of hadron h, and denotes its "cr"-term (Yih/m q may 
be interpreted as the number of qq pairs inside hadron h which diminish the (negative) 
qq density of the condensate). In obtaining Eq. ([8]), the GOR relation ([5]) has been used. 
Alternatively, one can directly use the definition of the quark condensate in terms of the 
quark-mass derivative of the thermodynamic potential, 

and evaluate the temperature and density-dependent part, Ct(fiB,T) = Q(fiB,T) — f2 V ac, in 
the free gas approximation. A similar strategy can be adopted for the gluon condensate, 
by utilizing its relation, Eq. (j3|), to the trace anomaly, 

((G 2 )) = -(e-3P)+m q ((qq)) , (10) 

and estimating the qb- and T-dependent parts of pressure and energy density in suitable 
expansions. At finite temperature, for a massless pion gas, one has e=3P and thus no 
correction to order T 4 (the system is scale invariant). It turns out that the lowest-order 
interaction contribution from (soft) tttt scattering does not contribute either so that the 
leading temperature dependence of the gluon condensate arises at order T 8 [53]. With the 
leading nuclear-density dependence as worked out in Refs. [55j[56], one has 

«G 2 » = -(G 2 ) - (m N - E N ) eN - |^ (ln^ - (11) 

(A p ~275MeV is a renormalization scale). The above relations allow for some interesting 
insights. As already noted in Ref. |56j . the linear-density expansions suggest that the gluon 
condensate is much less affected then the quark condensate, cf. also the upper panels in 
Fig. |3J It is not obvious whether recent finite-T lattice computations in QCD with 2+1 
flavors support this picture (see the lower panels of Fig. [4|), especially when approaching 
the critical temperature: both ((qq)) and ((G 2 )) drop significantly around T c and reach 
approximately zero at roughly 1.5 T c (the perturbative interaction contribution to e — 3P 
renders the gluon condensate negative at high T) . The low-density expansion of the quark 
condensate seems to suggest that temperature effects are weaker than density effects (upper 
left panel in Fig. |3J. This is, in fact, not the case: as a function of (pion-) density, the 
leading reduction of ((qq)) in a heat bath is quite comparable to cold nuclear matter, as 
determined by the coefficient in Eq. ([8]) which is in essence given by the respective a terms, 
£^70 MeV compared to £tv=45±15 MeV. 

Another interesting observation can be made when taking the expectation value of the 
trace of the energy momentum tensor over a single nucleon state, 

(N\T£\N) = -m N = (N\G 2 \N) + m g (N\qq\N) . (12) 

Since the second term (related to the a term) is small (or zero in the chiral limit), this 
relation seems to suggest that the major part of the nucleon mass is generated by the 
gluon condensate. This is to be contrasted with effective quark models (e.g., Nambu 
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Figure 4: Upper panels: density and temperature dependence of the chiral (left) and the 
gluon condensate (right), normalized to their vacuum values \57$ , obtained in a low-T and - 
qn expansions. Lower panels: T dependence at [i q =0 of the (subtracted and renormalized) 
chiral condensate (left) \5E^ and of the gluon condensate (right) f5P) / as obtained from 
Nf=2+1 lattice QCD computations. 



Jona-Lasinio) which, in mean field approximation, attribute the constituent quark mass 
entirely to the quark condensate, m*=G e s{qq). One should also note that, at least in 
the QGP phase, the vanishing of the quark or gluon condensate does not necessarily 
imply quark "masses" to vanish. E.g., in perturbative QCD, partons in the QGP acquire 
a thermal mass rrigg^gT. This mass term does not break chiral symmetry (its Dirac 
structure includes a 70 matrix) and presumably persists until close to T c , thus supplanting 
the constituent quark mass, m*, well before the latter vanishes. This has, of course, 
important consequences for the masses of hadronic states in the vicinity of T c . In addition, 
large binding-energy effects can be present, e.g., for the pion: if the chiral transition is 
continuous, the pion's Goldstone-boson nature could very well imply that it survives as a 
bound state at temperatures above T c [60l WD, [62] . 

As emphasized above, the only known direct way to extract observable consequences 
of changes in the QCD condensate structure is to probe medium modifications in its 
excitation spectrum. This applies in particular for the quark condensate which has a 
rather small impact on the bulk properties of QCD matter being suppressed by m q (the 
relation of the gluon condensate to the equation of state could, in principle, be tested via 
hydrodynamic or transport properties, but this turned out to be difficult in the context 
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of heavy- ion collisions [ID])- The generic model-independent consequence of xSR for the 
in-medium hadronic spectrum is the degeneracy of the spectral functions within chiral 
multiplets (i.e., for chiral partners) , e.g., vr-V, iV-A*(1535) and p(770)-cti(1260). In the 
"<t" channel (which asymptotically corresponds to a scalar- isoscalar pion pair), interesting 
medium effects have been observed in pion- and photon-induced production of S-wave pion 
pairs off nuclei [63\ [641 [65] . An accumulation of strength close to the two-pion threshold 
(which is not observed in the isotensor tt-tt channel) has been associated with an in- 
medium reduction of the "a" -meson mass as a precursor effect of xSR [66] (note that the 
leading-density approximation, Eq. ([8]), predicts a reduction of the quark condensate by 
~30% already at normal nuclear matter density). However, nuclear many-body effects |67[ 
168] . in particular the renormalization of the pion propagator in the nuclear medium, can 
essentially explain the experimental findings This raises an important question: to what 
extent do "conventional" in-medium effects encode mechanisms of xSR? From the point of 
view of the "<r" spectral function alone, it is not possible to distinguish whether a softening 
is caused by many-body effects or genuine mass changes figuring via medium modifications 
of the mass parameter in the underlying effective Lagrangian. Thus, a distinction of 
medium effects into "conventional" ones and those associated with an apparent "direct" 
connection to x$R is meaningless. Rather, a careful and exhaustive treatment of hadronic 
many-body effects is an inevitable ingredient for evaluating mechanisms of xSR- As already 
alluded to in the Introduction, a practical problem of using the tttt decay channel for 
studying medium effects are the strong final-state interactions of the individual pions when 
exiting the nuclear medium [69]. The same applies to the heavy-ion collision environment, 
implying that the tttt channel can only probe the dilute stages of the produced medium. 
This problem is overcome by dilepton final states, on which we will focus in the following. 



3 Vector Mesons in Medium 

3.1 Dileptons and Electromagnetic Correlation Function 

For a strongly interacting medium in thermal equilibrium the production rate of dileptons 
can be cast into the form [16] [T8], 

^ = f B ^;T) \ Blv lmUZ(M, q ^ B ,T) . (13) 

This expression is to leading order in the electromagnetic (EM) coupling constant, a ertl , 
but exact in the strong interaction. The latter is encoded in the EM spectral function, 
defined via the retarded correlator of the hadronic EM current, Jem (2^), 

KZ(%,q)) = -iJ d 4 x tP* e(z ) «^(x),jSa(0)]» • (14) 

In the vacuum, the spectral strength is directly accessible via the total cross section for 
e + e~ annihilation, 

a{e + e~ hadrons) = -^2i ( ~ 12?r) I m n™ c (s) , (15) 

s s 

recall Fig. [T] (the first factor is simply a(e + e~— >fi + fi~)=4Tra 2 m /s). As a function of invari- 
ant dilepton mass, M 2 =q L — q 2 , the spectrum basically decomposes into two regimes. In 



4 Similar results are obtained from a transport treatment of pion reinteractions in the medium [69] , 
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the low-mass region (LMR, M<1 GeV), the strength is absorbed in the three vector mesons 
p(770), uj(782) and 0(1020) representing the lowest resonances in the two-pion, three-pion 
and kaon-antikaon channels, respectively. Thus, the EM current is well described within 
the vector dominance model (VDM) [19] as given by the field current identity, 



m, 



m: 



nr 



j£ m (M < 1 GeV) = -V + + -£<, 



9 P 



9u 



(16) 



In the intermediate mass region (IMR, 1 GeV<M<3 GeV), the strength is reasonably well 
accounted for by a partonic description, 



j» m (M > 1.5 GeV) 



q=u,d,s 



e„ q^q 



(17) 



where e q denotes the electric quark charge in units of the electron charge, e. The connection 
between the two representations can be exhibited by rearranging the charge-flavor content 
of the quark basis into hadronic isospin quantum numbers, 



Jem = ~j= V>7^ 



uu — dd 1 uu + dd \/2 _ 
~7T + ?>~J2 3~ SS 



(18) 



reflecting the quark content of the (normalized) p (isospin 1=1), u) (1=0) and </> (1=0) 
wave functions, respectively. Converting the isospin coefficients into numerical weights in 
the EM spectral function, one obtains 



im n„ 



1 2 
lmD fJ + - Im D w + - Im Da 
9 9 



(19) 



which identifies the isovector (p) channel as the dominant source (experimentally it is 



even larg er as given by the electromagnetic decay widths, Tp — >ee/^uj — > 
evaluating the EM correlators using the currents (fT6|) and (fT7|) yields 



11). Explicitly 



Imn- C (M) 



Im£>^ ac (M) 



/ 1 _i_ <*s{M) , 
12-K 7T _r 



Nr. 



M <r A/f vac 



(20) 



q=u,d,s 



dual 



1.5 GeV, N c =3: number of quark colors, Dy = 1/[M Z 



■Sv] : vector-meson 



propagators). The associated processes in the thermal dilepton production rates are, of 
course, the inverse of e + e~ annihilation, i.e., two-pion, three-pion and KK annihilation 
(channeled through the p, uj and (j)) in a hadronic phase! and qq annihilation in a QGP. 
But what about hadronic emission in the IMR and QGP emission in the LMR? The 
former follows from time-reversal invariance of strong interactions: to the extent that the 
hadronic final state in e + e~ annihilation can be represented by a statistical (thermal) 
distribution (which is empirically approximately satisfied), hadron-gas emission in the 
IMR corresponds to multi-hadron annihilation (4-7T, 6tt — > e + e~, etc., which may be built 
from 2/9, nai, iroj, etc.), with a total strength given by the partonic continuum. QGP 
emission in the LMR is, of course, closely related to a central question of this review: How 



5 Note that the dominance of the isovector channel is naturally associated with the annihilation of the 
two lightest constituents in a hadronic medium. 
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Figure 5: Left panel: vector and axialvector spectral functions as measured in hadronic 
t decays J57J/ with model fits using vacuum p and a± strength functions supplemented 
by perturbative continua |70| /; right panel: scenarios for the effects of chiral symmetry 
restoration on the in-medium vector- and axial-vector spectral functions. 



does the dilute hadronic resonance gas rate evolve into the chirally restored, deconfined 
QGP rate? At sufficiently low temperatures and/or baryon densities virial expansions in 
a hadronic basis can provide initial insights. With increasing T and qb resummations 
become necessary for which many-body approaches are a suitable tool. It is currently 
an open question how far up in qb and T these calculations are reliable. Selfconsistent 
schemes are, in principle, capable of describing phase-transition dynamics, which, ideally, 
could be constrained by unquenched lattice-QCD calculations of the dilepton rate below 
T c (energy sum rules turn out to be particularly useful to connect spectral functions to 
order parameters). Eventually, in the high-temperature limit, the LMR rate should recover 
perturbative qq annihilation, where a systematic evaluation of corrections becomes feasible 
again. The remainder of this section is devoted to a discussion of these approaches. With 
hindsight to Sec. [H we will focus on models for which quantitative applications to dilepton 
observables have been made, with the isovector {p) channel playing the leading role. In the 
following, for brevity, we refer to the IJ P =11 ± chiral partner channels as vector (V) and 
axialvector (A) ones. In the vacuum, both can be well represented by a low-lying resonance 
pole (p and a±) and a continuum above, see left panel of Fig. Two schematic scenarios 
for the degeneration of vector and axialvector channels at chiral restoration ("dropping 
mass" and "resonance melting" ) are sketched in the right panel of Fig. [5j 

3.2 Medium Effects I: Model Independent 

In principle, model-independent assessments of medium effects do not involve free param- 
eters. These can be realized by virial expansion schemes based on experimental input for 
vacuum spectral functions (valid for dilute hadronic matter), perturbative QCD calcula- 
tions (valid in the high-T limit) or first-principle lattice QCD computations. 
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3.2.1 Chiral Reduction and Mixing 

The leading temperature dependence of vector and axialvector correlators, i.e., their 

modification in a dilute pion gas, can be inferred from chiral reduction and current algebra. 
They allow to simplify 1-pion matrix elements of any operator according to 

(01010 = "f (0\m,O\\0) , [Q%j^ A ]=^ abc fAy , (21) 

Jtt 

where {a, b, c} are isospin indices. Evaluating the Fourier transforms of the thermal expec- 
tation values in the chiral and soft pion limit (i.e., m n =0 and neglecting any momentum 
transfer k from thermal pions in the heat bath), one obtains the "mixing" theorem |71j 

U v>A (q) = (1 - e) n^(g) + e n^(g) (22) 

with the mixing parameter e=T 2 /6f% (the Lorentz structure remains as in the vacuum). 
The leading-T effect on the V and A correlators is a mere admixture of the chiral partner 
with a corresponding reduction of its original strength, via processes of the type ir+V <-> A 
and 7r + A <-> V; width and mass of the vacuum correlators are unaffected. For dilepton 
production, this implies a reduced p pole strength as well as an enhancement of the "dip" 
region, M~l-1.5GeV, where the a\ resonance provides a "maximal feeding". 

When naively extrapolating the mixing expression, Eq. (j22[) . to chiral restoration 
(e:=l/2), one finds T c =\/3/7r=160 MeV. This is, however, misleading for several reasons. 
First, this estimate does not coincide with a similar extrapolation for the vanishing of the 
chiral condensate, cf. Eq. (|8|). Second, even a moderate amendment in terms of a finite 
pion mass in the scalar density shifts the estimate to T c ~225 MeV. Both facts underline the 
inadequacy of the extrapolation of a lowest-order result. Third, the chiral and soft-pion 
limits are kinematically not a good approximation (e.g., at T=150MeV, thermal pions 
typically bring in an energy of ~300-400 MeV). In cold nuclear matter, a similar mixing 
is operative via the coupling of the pion cloud of p and a\ to the nuclear medium [72 \ 173]. 

A much more elaborate treatment of the chiral reduction formalism has been conducted 
in Refs. [7H [T5j [76] . These calculations are based on realistic fits to vacuum correlators, 
do not invoke kinematic approximations (chiral or soft-pion limits) and include both pion 
and nucleon ensembles. The leading-density part has been subjected to constraints from 
nuclear photo-absorption including the first and second resonance region via A(1232) and 
A^(1520) excitations. This allows for meaningful applications to dilepton spectra which 
have been carried out and will be discussed in Sec. [H Note that these calculations do 
not explicitly invoke the notion of VDM, but the fact that the vacuum correlators are 
constructed with p and a\ pole dominance, which is not upset in the linear density scheme, 
implies that VDM is still present upon inclusion of medium effects. 

3.2.2 Lattice QCD and Susceptibilities 

First-principle computations of light-hadron correlation functions in medium are based 
on a lattice discretized form of the finite-T QCD partition function. Besides a finite lat- 
tice spacing, additional approximations currently involve the restriction to finite volumes 
as well as the use of unphysically large up- and down-quark masses in the simulations. 
Furthermore, the implementation of chiral symmetry is not trivial in the lattice formula- 
tion. The numerical evaluation of the QCD path integral, is facilitated by transforming 
the action to imaginary (Euclidean) time, which converts the oscillatory behavior of the 
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Vector 




Figure 6: Left panel: vector correlation function as a function of Euclidean time as eval- 
uated in quenched lattice QCD for a gluon plasma at temperatures above T c J77| /. The 
in-medium correlators are normalized to the free one using the integration Kernel at the 
same temperature (the so-called "reconstructed" correlator). Right panel: thermal dilepton 
rates, dN/(d 4 qd 4 x), in quenched lattice QCD as extracted from the correlation functions 
shown in the left panel using the maximum entropy method. The lattice results are com- 
pared to calculations in perturbation theory, either to leading order (0(a®) ) qq annihilation 
(solid line) or within the hard-thermal-loop (HTL) framework J7£| / (dashed line). All rates 
are calculated at a total pair 3-momentum of q=0, i.e., the dilepton energy, co=qo, equals 
its invariant mass, M . 



integrand in the partition function into an exponential damping. The pertinent Euclidean 
correlation function, Il(r), is related to the physical spectral function, p = — 2Imil, via 

oo 

tt( t\ f d( & < t\ cosh[( g0 (r-l/2T)] 

n(T ' J ~2tt ' ?! T) sinh[W2T] " (23) 
o 

The resulting Euclidean vector correlators in "quenched" above T c shows a moderate 

enhancement over the free correlator, cf. left panel of Fig. [6] [77J. The extraction of 
the spectral function requires an inverse integral transform over a finite number of r 
pointaj which can only be achieved with a probabilistic treatment based on the "Maximum 
Entropy Method" [79] . The resulting strength function has been inserted into the dilepton 
rate and is compared to perturbative QCD (pQCD) rates in the right panel of Fig. [6l The 
leading-order pQCD corresponds to the qq strength distribution in Eq. ([20]) . lower line, 
while the hard-thermal-loop (HTL) improved rate is from Ref. [78]. The latter shows the 
expected divergence for M — ► which is caused by the Bose factor and photon propagator 
which overcome the p oc q§ dependence of a retarded correlation function (cf. also Ref. [80]). 
This feature is not shared by the lattice result which might be an artifact of, e.g., the finite- 
volume restriction (it would also suggest a small or even vanishing photon production 
rate). On the other hand, the enhancement in the Euclidean correlator translates into an 
enhanced dilepton rate at energies of a few times the temperature. Whether this reflects 
a broad resonance structure is not clear at present. 



6 In the "quenched" approximation the fermionic part of the QCD action is neglected in the evaluation 
of the Euclidean path integral. This amounts to neglecting fermion loops. 

7 The (anti-) periodicity of the boson (fermion) fields at finite T restricts the Euclidean time direction 
to the interval [0, 0\ where (3 = 1/T. 
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Figure 7: Isoscalar (left) and isovector (right) quark-number susceptibility for various 
quark chemical potentials, p, q = pb/3, as computed in unquenched lattice QCD 181)/ . 

Additional constraints from lattice QCD are provided by susceptibilities which are 
defined as second-order derivatives of the thermodynamic potential. In our context, the 
quark- number susceptibilities are of special interest, 

Xa~ jrY-ILaiqo = 0,(7^0) , (24) 

which can be decomposed in isoscalar (p q = (p u + pd)/2) and isovector (pj = (p u — 
Pd)/2) channels carrying the quantum numbers of the to and p, respectively. The spacelike 
limits of the correlators basically represent the screening masses in the respective channels. 
Lattice QCD computations of the quark-number susceptibilities indicate that both p and 
<jj channels behave smoothly with temperature for small chemical potentials, see Fig. [7J 
However, as p q increases, Xq develops a peak whereas xi remains smooth. The former 
indicates an increase in the (local) baryon-number fluctuations and may be a precursor of 
the baryon-number discontinuity between hadronic and QGP phase as one is approaching 
a first-order line. Remarkably, this is not seen for the isospin fluctuations. 

3.3 Sum Rules and Order Parameters 

Sum rules are currently the most promising tool to connect the nonperturbative physics 
encoded in spectral functions to the condensate structure of the QCD vacuum. In partic- 
ular, the Weinberg sum rules directly relate order parameters of xSR to the axial-/vector 
spectral functions, which, in the medium, have not been exploited much to date. 

3.3.1 Chiral Sum Rules 

The Weinberg and DMO sum rules [52|. [53] . Eqs. ([6]), directly relate moments of the 
"vector minus axialvector" spectral functions to chiral order parameters. This is a rather 
fortunate situation in view of the dominant role that the isovector- vector (p) channel plays 
in dilepton production, recall Eq. (|19p . For Nf=2, the ui is a chiral singlet, while in the 
strangeness sector ((ft), i.e., for Nf=3, chiral symmetry becomes much less accurate (e.g., 
{{ss)) persists much farther into the QGP). 
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As has been shown in Ref. [82J, the Weinberg sum rules remain valid at finite tem- 
perature, albeit with two important modifications induced by the breaking of Lorentz 
invariance caused by the heat bath which defines a preferred rest frame: (i) each energy 
sum rule applies for a fixed three-momentum, and (ii) at finite three-momentum, the 
vector and axialvector spectral functions split into longitudinal and transverse modes, 

= n^p,f + n^pf , (25) 

with individual sum rules for each of them. The explicit form is as follows: 

00 

= 0, (26) 

= 0, (27) 

= -27ra,((0 4 )) • (28) 
n 

In writing Eqs. ()26|) - (j28[) the pionic piece of the (longitudinal) axialvector correlator has 
been absorbed into the definition of the in-medium spectral function, ImlI^(go, q)', in the 
vacuum and in the chiral limit it is represented by a sharp state, lmIL% l '=f 2 M S(M )P£ V . 
In this form it only contributes to the first sum rule, Eq. (|26f) . However, in matter (and 
for m 7r >0) this is no longer true since the pion is expected to undergo substantial medium 
effects. 

The in-medium chiral sum rules constitute a rich source of constraints on both energy 
and three-momentum dependence of in-medium spectral functions. The energy moments 
demonstrate that chiral restoration requires degeneracy of the entire spectral functions. 
Combining 1QCD computations of order parameters with effective model calculations thus 
provides a promising synergy for deducing chiral restoration from experiment [43j . 



/ < t % 2^ [Imll£(?0,9) " lmU^(q ,q) 
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3.3.2 QCD Sum Rules 

QCD sum rules are based on a (subtracted) dispersion relation for a correlation function 
in a given hadronic channel a, formulated for spacelike momenta q 2 =—Q 2 <0 [83], 

n Q (Q 2 ) = n Q (o) + w a (o) q 2 + q 4 / . (29) 

J irs z s + Q z 

The right-hand-side (rhs) contains the spectral function which is usually related to ob- 
servables or evaluated in model calculations. On the left-hand-side (Ihs), the correlation 
function is expanded into a power series of 1/Q 2 (operator-product expansion = OPE) 
where the (Wilson) coefficients contain perturbative contributions as well as vacuum- 
expectation values of quark and gluon operators (the nonperturbative condensates; for 
practical purposes the convergence of the OPE is improved by means of a so-called Borel 
transformation which we do not discuss here). The explicit form of the OPE for vector 
and axialvector correlators reads (truncating higher order terms in m q , a s , etc.) 

Q 2 ~ Svr 2 ln v 2 + Q* + 24Q4 7r ^ v 81Q6 W> + ^ 

Q 2 ~ 87T 2 ^ Q4 + 24Q 4 8 1Q6 W + ^ 
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Figure 8: QCD sum rule constraints on p meson mass and width as inferred from Breit- 
Wigner parameterizations of its spectral function Jff7| / (left panel: vacuum, right panel: cold 
nuclear matter at saturation density). "Allowed" regions of mass and width are indicated 
by the bands between solid and dashed curves, corresponding to maximal deviations between 
the l.h.s and r.h.s. of the SR of 0.2% and 1%, respectively. The diamond depicts the 
corresponding vacuum parameters. 



where the four-quark condensates have been approximated by factorizing them into the 
squared two-quark condensate with parameters k, k which simulate intermediate states 
other than the ground state (the scale p is typically chosen around 1 GeV). Note that SBxS 
is nicely reflected by the opposite signs of the quark-condensate terms in II y and II^ , while 
the "flavor-blind" gluon condensate enters with the same sign. Qualitatively, the (positive) 
gluon condensate actually induces a softening of the spectral function (i.e., a larger weight 
at small s in the dispersion integral) [84j. On the other hand, for the vector channel, the 
negative contributions from the quark condensates on the Ihs of the sum rule push spectral- 
function strength to larger s, relative to the axialvector channel (this may seem surprising 
in view of the masses of the pertinent resonances, m p =0.77GeV vs. m ai =1.23 GeV; recall, 
however, that the (longitudinal) axialvector channel contains a contribution from the 
axialvector current of the pion). Inserting numerical values, a s =0.35, m g =0.005 GeV, 
(w)=(-0.25GeV) 3 and (^G^ 2 )=0.012 GeV 4 , leads to 

ny? c If O 2 0.0062 GeV 4 0.039 GeV 4 0.029 GeV 6 k\ 

■l.llln^ — + — , (32) 



Q 2 8tt 2 V ' m 2 Q 4 Q 4 Q 6 

illustrating that the leading contributions arise from the gluon and four-quark condensates 
(especially for typical values of k~2.5), while the impact of the quark condensate is rather 
moderate. For the vector channel, and in vacuum, there is a large cancellation between 
the gluon and 4-quark condensate terms. However, in the medium this is presumably 
lifted, especially at low T and pb where quark and gluon condensates change rather 
differently. The stronger reduction of the "repulsive" 4-quark condensate relative to the 
"attractive" gluon condensate induces a softening of the spectral function in the dispersion 
integral. The softening can be satisfied by both broadening and/or a downward mass 
shift [85l EH [87], EEl [89]. Quantitative studies (which also include effects of non-scalar 
condensates induced by hadron structure of the heat-bath particles) based on Breit-Wigner 
model spectral functions are displayed in Fig. [8l For the axialvector channel, the reduction 
in both condensates suggests a substantial loss of soft-mode strength which points at the 
dissolution of the pion mode (whose polestrength is given by /„-) as a consequence of (the 
approach toward) xSR- 
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Finally, it is instructive to compare uj and p mesons: while their OPE side is rather 
similar (governed by 4-quark condensates), the subtraction constant, Ily(O) = qn/^Mn 
to leading order in qn, makes a difference. It is given by the Thompson limit of the VN 
scattering amplitude and turns out to be identical in the p and u sum rule. However, 
since Im Yl p is larger than Im n w by an isospin factor of (g^/gp) 2 — 9 (recall Eq. (|19p ). the 
finite- qn subtraction actually stabilizes the uj sum rule, implying stronger medium effects 
(softening) on the p than on the uj (it amounts to a "repulsive" contribution on the OPE 
side counterbalancing the reduction in the 4-quark condensate). 

3.4 Medium Effects II: Chiral Effective Models 

Model- independent and/or low-density approaches as discussed above provide valuable 
constraints on the vector and axialvector correlators and their connections to QCD vac- 
uum structure. However, quantitative calculations suitable for comparison with experi- 
ment require the construction of effective models. As indicated in the Introduction, in 
the low-mass region most of the thermal dilepton yield in heavy-ion collisions is expected 
to emanate from the hot / dense hadronic phase (even at collider energies) , especially from 
the p channel. Hadronic chiral Lagrangians are therefore a suitable starting point, ex- 
tended by the implementation of the low-lying vector mesons. This is usually done by a 
local gauging procedure of the chiral pion Lagrangian, thus realizing the gauge principle 
at the composite (hadronic) level. The most common approaches are based on non-linear 
realizations of chiral symmetry (i.e., without explicit a meson) within the Hidden Local 
Symmetry (HLS) [90] or Massive Yang Mills (MYM) [91] schemes. Rather than review- 
ing these in a comprehensive form, we here focus on recent developments with relevance 
for dilepton production, i.e., the "vector manifestation" (VM) scenario of SB^S within 
HLS [92] fSec. [3XT]) . as well as hadronic many-body theory within MYM (Sec. [3X2} . 

3.4.1 Hidden Local Symmetry and Vector Manifestation 

Within the HLS framework, an alternative realization of chiral symmetry in the meson 
spectrum has been suggested in Ref. |92| . by identifying the chiral partner of the pion 
with the (longitudinal) p (rather than with the a). This "vector manifestation" of chiral 
symmetry has been shown to give a satisfactory phenomenology of hadronic and EM decay 
branchings in the vacuum. When applied within a finite-T loop expansion, the p-meson 
mass was found to be affected at order T 4 (consistent with chiral symmetry), showing a 
slightly repulsive shift. However, when matching the hadronic axial-/vector correlators to 
pQCD in the spacelike regime (using an OPE), a reduction of the bare p mass has been 
inferred, consistent with "Brown- Rho" scaling [21]. In addition, vector dominance was 
found to be violated in the medium, leading to a gradual decoupling of the p from the 
EM current toward the critical temperature. However, the finite-T EM formfactor, which 
determines the dilepton production rate, clearly shows the downward moving p peak [93J, 
see Fig. [9j An interesting question is how these features develop in the presence of finite 
baryon density. 
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Figure 9: Pion EM formfactor at finite temperature in the HLS vector manifestation 
framework 193)/ : left panel: with finite-T loop effects; right panel: additionally including a 
T -dependence of the bare p mass above T=0.7T C . 



3.4.2 Massive Yang-Mills and Hadronic Many-Body Theory 

As in HLS, the basic building block of the MYM Lagrangian is the chiral pion Lagrangian 
based on the unitary pion field, 



a V2 



U = exp(fv / 2<A// 7 r) , 
Hadronic gauge fields, A^ L R are introduced via the covariant derivative, 



(33) 



D»U = ig(A£U - UA R ) (34) 
and supplemented with kinetic and mass terms (with bare mass mo). One has 

C, 



■-mym 



\fl tr \D»UD»UA - \ tr [(Ff ) 2 + (Fj? ) 2 ] + m 2 , tr [(A»f + «) 2 ] 



-i£ tr 



+ a tr 



F^UFr^uA , (35) 



where the last two (non-minimal) terms are necessary to achieve a satisfactory phe- 
nomenology in the vacuum. After the identifications = V^=A R + A^, A fl =A R — A^ L 
(and a field redefinition of the axialvector field to remove a d^nA^ term), the leading 
terms of the MYM Lagrangian take the form 
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mym 
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2 r -2 -2 
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(36) 



(gp n7T =^g ). Note that the Higgs mechanism induces the splitting of p and a\ masses, 
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which is entirely due to SB^S (via /„■). The bare p mass itself is an external parameter, 
which is different from the HLS scheme discussed in the previous section. Electromag- 
netism is readily included into the MYM Lagrangian by adding the vector dominance 
coupling [91] 
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Figure 10: Graphical representation of self-energy diagrams characterizing the interactions 
of the p meson in hot and dense hadronic matter: (a) renormalization of its pion cloud due 
to modified pion propagators, and direct interactions of the p meson with (b) baryons and 
(c) mesons, typically approximated by baryon- and meson-resonance excitations \94[ \3C 



where B p denotes the photon field. In this scheme, VDM remains valid in the medium, 
and the task of computing the low-mass isovector axial/-vector correlators amounts to 
assessing the medium modifications of p and a\ mesons. 

The />meson propagator in hot and dense hadronic matter can be written as 

D L p ' T (q , q; HB,T) = * yL,T y L,T > ( 39 ) 

with transverse and longitudinal modes as defined in Eq. (|25p . The key quantities are 
the in-medium selfenergies, Xp' T , which may be classified as follows (cf. Fig. I1U|) : Sp T7r 
accounts for the pion cloud of the p, which in the vacuum gives rise to its finite width 
via p — ► 7T7T. Direct interactions of the p with mesons (M=ir, K, p, . . . ) and baryons 
(B=N, A, A, ... ) from the heat bath are represented by £ p m and S p _b, respectively; they 
vanish in the vacuum. In terms of underlying scattering processes, the latter are typically 
resonance excitations (e.g., pir — > a\ or pN — > A r (1520)) while medium modifications 
of pions (e.g., ttN — > A) in S p7r7r correspond to, e.g., t-channel ir exchange processes 
(pN — > 7rA). When evaluating interactions which are not directly constrained by chiral 
(or gauge) symmetry (especially those involving higher resonances), phenomenological 
information is essential for a reliable determination of the coupling constants (and cutoff 
parameters in the hadronic formfactors to account for the finite size of the hadrons). The 
simplest form of such constraints are hadronic decay widths of resonances (e.g., a\ — > pir), 
supplemented by radiative decays (e.g., a\ — > / yir). However, especially for "subthreshold" 
states (e.g., uj — > pir or iV(1520) — > pN), where the coupling is realized via the low-energy 
(tttt decay) tail of the p spectral function, empirical information can be rather uncertain. 
In this case, comprehensive constraints inferred from scattering data become invaluable. 
Unfortunately, in practice this is only possible for pN interactions (e.g., via irN — > pN or 
jN scattering), but, as it turns out, the modifications of the p due to interactions with 
nucleons are generally stronger than with pions. In addition, by using nuclear targets, one 
has the possibility to constrain (or test) the modifications in nuclear matter, rather than 
on a single nucleon (which corresponds to the leading-order density effect). 

Let us start by discussing finite-T effects. Calculations of the p propagator in a 
hot pion gas based on the MYM scheme [95] have shown small medium effects. An 
extended analysis [96] of the p in hot meson matter, including resonance excitations 
(pir — > ai, uj, h\, it', a2(1320), pK — > K*,K\, pp — > /i(1285)) and pion Bose enhance- 
ment in S p7r7r , leads to total broadening of ~80MeV at T=150MeV (corresponding to a 
pion density ^=0.12 fm~ 3 ~0.75 Qo), with little mass shift. Approximately ~20MeV of 
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q 2 [GeV 2 ] 

Figure 11: Comparison of p-meson spectral functions in cold nuclear matter within the 
hadronic many-body approaches of Refs. 1100)/ (left panels, based on two different phase 
shift analysis of ttN scattering }101\ 102]) and \94\ \ 30 $ (right panel). 

the broadening is due to the tttt Bose factor (cf. also Ref. [97]) and ~50-60MeV due to 
meson resonances. The latter is comparable to Refs. \98\ [99] which are directly based on 
pir and pK scattering amplitudes. 

Next, we turn to modifications in cold nuclear matter. Fig. QT] shows two calcula- 
tions in which the underlying p self-energies have been rather thoroughly constrained. 
In Ref. [100] (left panels), a pN resonance model (corresponding to £ p at) has been con- 
structed utilizing a detailed analysis of empirical ttN — ► pN phase shifts and inelastici- 
ties |10H 1102] , The resulting p spectral functions are displayed at normal nuclear density 
(taken as £>tv=0.15 fm -3 ) for various three-momenta and two distinct data sets for con- 
straints. A substantial broadening of close to ~200MeV is found, with a slight upward 
peak shift of a few tens of MeV; the three-momentum dependence is relatively weak. In 
Refs. [Ml GO], S p 7v and an in-medium pion cloud, S p7r7r (incorporating P-wave "pisobar" 
nucleon- and A-hole excitations and associated vertex corrections), have been calculated 




M (GeV) M (GeV) 

Figure 12: Comparison of p-meson spectral functions in hot and dense hadronic mat- 
ter from Refs. (dashed lines) and Wty (solid lines) at temperatures of T =100 MeV 
(left panel) and T =150 MeV (right panel) for nucleon densities of qn=0.5,1.0 Qo (the 
corresponding nucleon chemical potentials are [in =673,745 MeV for T =100 MeV and 
fi N =436,542 MeV for T =150 MeV, respectively). 
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Figure 13: In-medium spectral functions of light vector mesons in hot and dense hadronic 
matter based on the hadronic many-body approach of Refs. \103[ [3ff\ 1 1 05^ : upper panel: 
under conditions resembling heavy-ion collisions at SPS (i.e., along an isentropic tra- 
jectory in the phase diagram which preserves the measured hadron ratios determined at 
(^f em ,r chem ) ~ (230,175) MeV) fTUbJ : lower panels: at (/x B ,T) ~ (25,180) MeV (re- 
sembling chemical freeze- out at RHIC) with (short-dashed lines) and without (long-dashed 
lines) medium effects induced by anti- /baryons 1201/ . 



and constrained by total photoabsorption data on the nucleon and nuclei |103j . as well as 
total irN — > pN cross sections. The resulting spectral functions are quite similar to the 
ones of Ref. [100] . with a somewhat stronger broadening of ~300MeV at ^=0.16 fm -3 
and a comparable mass shift of ~40MeV. It is quite remarkable that the predicted in- 
medium mass and width of ~ (810,450) MeV are in good agreement with the QCD sum 
rule constraints derived in Ref. [87], cf. right panel in Fig. [8] Both broadening and mass 
shift decrease at higher three-momentum, e.g., (Am p , Ar p )~(30,150) MeV at g=lGeV. 
Both calculations [1001 [94] include a rather strong coupling to /?iV(1520) iV™ 1 excitations 
(appearing as a low- mass peak or shoulder in the p spectral function). This has been ques- 
tioned in Ref. [104] based on a coupled channel analysis of S'-wave pN and toN scattering, 
where all nucleon resonances but the A(1232) are generated dynamically via four-point 
interactions. The (generated) 7V(1520) is deduced to primarily couple to ojN rather than 
pN, entailing an in- medium p with significantly less broadening. 

Finally, we turn to a hot and dense hadronic medium as expected to be formed in high- 
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energy heavy-ion collisions. In Fig. [T2] the p spectral functions of the hadronic many-body 
calculations [944 [30] are compared to those obtained in the scattering-amplitude approach 
of Ref . . The latter exhibit less broadening and a small (if any) downward mass shift of 
the p peak, compared to the upward shift in the many-body approach (mostly induced by 
baryonic effects). Consequently, in terms of spectral strength, the discrepancies between 
the two calculations are largest for masses around M~0.7GeV, as well as for very low 
mass, M<0.4 GeV, by up to a factor of ~2. This mass region is much magnified in thermal 
dilepton production rates due to the Boltzmann factor and a photon propagator oc 1/M 2 . 
However, the amplitude approach only accounts for interactions with pions and nucleons, 
while the many-body calculations include estimates of p interactions with strange baryons 
and resonances [30]. This difference may account for some of the discrepancy. 

In preparation for applications to dilepton spectra in URHICs, we summarize in Fig. [13] 
in-medium l/-meson spectral functions in the many-body approach under conditions rel- 
evant for SPS (upper panels) [30 ] 1106] and RHIC (lower panels) [105] . The p meson (left 
panels) "melts" when extrapolated to temperatures close to the expected phase boundary. 
Baryons play an essential role in the melting, even at RHIC (where the net baryon density 
is small), since the relevant quantity is the sum of baryon and antibaryon densities. The 
effects due to baryons and antibaryons are most prominent as an enhancement in the mass 
region below M~0.5 GeV. The u> and especially 4> spectral functions (right panels) appear 
to be more robust. One should also point out that at T=120,150 MeV in the upper panels 
appreciable pion and kaon chemical potentials are present which sustain larger hadron 
densities and thus support stronger medium effects than in chemical equilibrium. 

3.5 Thermal Dilepton and Photon Rates 

The in-medium vector-meson propagators discussed in the preceding section are converted 
to thermal dilepton rates via Eqs. ([20]) (upper line) and ([T3]) . This is based on the as- 
sumption that VDM for the EM correlator remains valid in the medium^]. The resulting 
three-momentum integrated thermal dilepton rates are summarized in Fig. PT4l The left 
panel, which displays the isovector channel, reiterates that the p resonance signal dis- 
appears from the mass spectrum as one approaches the putative phase boundary. The 
hadronic rates also include an estimate of the leading-T chiral mixing effect, Eq. (122j) . in 
the mass region M=l-1.5GeV. The comparison to perturbative qq annihilation reveals 
that the the top-down extrapolated QGP rate closely coincide with the bottom-up ex- 
trapolated in-medium hadronic one, especially in case of the HTL-improved qq rate. This 
feature suggests that the hadronic rate has indeed approached xSR (since the QGP rates 
are chirally symmetric at any finite order in perturbation theory) [34[ [30] . The "matching" 
of QGP and hadronic rates occurs directly in the timelike regime without the need for 
in-medium changes of the bare parameters in the effective Lagrangian. Medium effects due 
to baryons play an important role in this mechanism; the situation is similar for small pb 
and close to T c where the sum of baryon and antibaryon densities is appreciable, see right 
panel of Fig. [TH uj and especially <j> mesons appear to be more robust, possibly surviving 
above T c . The dilepton rates in the vector manifestation of HLS [93] look rather different; 
based on the pertinent pion EM formfactor, Fig. [9] a distinct p peak survives in the rate 
up to temperatures of at least T=0.85 T c ~155 MeV (assuming T c ~180MeV). 

8 Strictly speaking, the EM correlator of Refs. [941130] includes corrections to VDM in the baryon sector 
as determined via photoabsorption spectra on the nucleon and nuclei 103 ; the assumption is that this 
modified version of VDM is not affected at higher densities and at finite temperature. 
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Figure 14: Three-momentum integrated thermal dilepton rates at fixed temperature |7Q| , 
\107l for the vacuum p (dotted lines), the hadronic many-body approach of Refs. \103[ 
(solid lines) and for the QGP using either free qq annihilation (dashed line) or hard- 
thermal loop improved rates |7&| / (dash-dotted line). The left panel refers to the isovec- 
tor (p) channel, under conditions resembling heavy-ion collisions at the SPS (fixed 
Pb=%70 MeV). The right panel additionally includes isoscalar (to and <p) channels and 
corresponds to small /j>b—%5 MeV appropriate for the conditions at collider energies. 



Emission rates of dileptons are closely related to those of real photons which are de- 
termined by the lightlike limit (go =| q\) of the EM spectral function, 

«°J^r = -^T f B (Qo;T) lmU cm (M = 0,q;p B ,T) . (40) 
a^xd^q 7r z 

In Ref. [108J the in-medium p of Refs. [103|, 130] has been found to constitute the dominant 
hadronic source of thermal photons for momenta up to g~l GeV; above, i-channel meson 
exchange reactions not included in the spectral function (most notably tt and lo exchange in 
irp — ► 7T7) take over, cf. left panel of Fig. [15j Similar to the dilepton case, at temperatures 
of 150-200 MeV, the strength of the combined thermal rate for hadronic photon production 
turns out be very comparable to perturbative QGP emission, especially for the complete 
leading-order result |110| . 

4 Interpretation of Dilepton Spectra 

In this section we will scrutinize experimental results for dilepton spectra in light of the 
theoretical developments elaborated above. A brief discussion of production experiments 
off nuclei, representing cold nuclear matter up to saturation density, will be followed by 
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Figure 15: Thermal photon production rates in hot hadronic matter 1108$ (left panel), and 
compared to 1-loop J 109$ and full leading-order 1110$ QGP emission (right panel). 
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a more extensive study of invariant-mass and momentum spectra in heavy-ion reactions 
involving hot and dense matter possibly probing the transition regime to the QGP. 

4.1 Medium Effects in Nuclei 

Dilepton production experiments off nuclei have the advantage over heavy-ion collisions 
that the medium is well-defined. Medium- to heavy ground-state nuclei resemble in their 
interior infinitely extended nuclear matter at vanishing temperature. Therefore the exper- 
iments probe to a large extent the properties of "cold nuclear matter". However, a good 
knowledge of the production process is required, and medium effects are typically rather 
moderate, further reduced by surface effects and decays outside the nucleus especially at 
large three-momenta (which, in turn, are needed in the production process). 

The E325 experiment at KEK [46] used 12 GeV proton projectiles and found significant 
differences in the spectra between C and Cu targets, see left panel of Fig. [TBI After 
subtraction of combinatorial background as well as r\ and u) Dalitz decays, the best fit 
to the excess spectra using co and p Breit-Wigner distributions was obtained with a mass 
shift of ca. 9% at nuclear matter density, and a p/to ratio of about ~0.45. 

The CLAS experiment [13] used a photon beam at Jefferson Lab with incident en- 
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ergies -E 7 =0.6-3.8 GeV, directed on various nuclear targets. After subtraction of the 
combinatorial background, determined with absolute normalization, the invariant-mass 
spectra are best reproduced with a p spectral distribution with a mass and width of 
(m p ,r p )=(775±5,220±15) MeV, cf. right panel of Fig. [T6j These values are well repro- 
duced by Boltzmann transport calculations and are consistent with the predictions of 
Refs. [Ml [30] at g N =V.b g and g=lGeV, where (Am p , Ar p )~(15,75) MeV, see Ref. [TT2] 
for a recent calculation. An apparent difference between the E325 and CLAS spectra is 
that the background subtraction in the former removes any excess for M~0.85-l GeV; this 
suppresses (and possibly shifts down) the p contribution in the E325 fit. 

Photoproduction experiments (i? 7 =0.8-1.12 GeV) of P-wave 7r + 7r _ pairs off 2 H, 3 He 
and 12 C were conducted by the TAGX collaboration |113j . The spectra for the 12 C target 
support medium effects in line with hadronic many-body p spectral functions [9U l30] . 



4.2 Heavy-Ion Collisions 

In contrast to production experiments off nuclei, the (energy-) density of the medium 
created in heavy-ion collisions undergoes a rapid evolution after initial nuclear impact 
until break-up. Even under the simplifying assumption of local thermal equilibrium, a 
good knowledge of the temperature and baryon-density evolution is necessary to convert 
the dilepton rates discussed above into a space-time integrated spectrum In addition, 
sources other than thermal radiation have to be considered, especially toward higher mass 
or qt where the assumption of equilibrium becomes increasingly questionable. These issues 
are addressed in Sec. 14.2.11 Phenomenological analyses of dilepton spectra focus on recent 
SPS data from NA60 and CERES/NA45 in Sees. K2?2\ and HT231 respectively, with a short 
digression to direct photons (W98) and a brief outlook to future experiments in Sec. 14.2.41 



4.2.1 Thermal Evolution and Dilepton Sources 

As discussed in Sec. II. 2[ hadronic observables in URHICs point at a reasonable degree 
of thermalization of the bulk medium produced in these reactions. Therefore, we here 
focus on hydrodynamic and expanding fireball approaches utilizing the assumption of 
local thermal equilibrium jf| 

Thermal emission spectra follow from the convolution of the production rate, Eq. fj 1 31) , 
over the space-time evolution of the medium, 

dN U M J, f l3 fd 3 q dN u 



I At [ d 3 x -j— |j— (M, q; T, p N , p s ,Pi) Acc{M,q t ,y), (41) 
J J i » d 4 x d 4 g 



dMdy Ay 

V FB 

where Acc(M,qt,y) accounts for the detector acceptance and Ay denotes the correspond- 
ing rapidity interval. The temperature and chemical potentials in general depend on 
space-time, (t,x). Note that while p^ and p s correspond to exact conservation of baryon 
number and strangeness, effective chemical potentials pi = P-k^k,^,... are needed to preserve 
the experimentally observed hadron ratios in the evolution of the hadronic phase between 
chemical (T ch =155-175 MeV) and kinetic freezeout (T fo =100-140 MeV). 

An overview of several key input parameters of three thermal approaches \36\ 11061 
11171 11181 1119^ 1120] . which have been used to compute dimuon spectra in comparison 



Comparisons of dilepton spectra computed in hydrodynamic/fireball and transport calculations (based 
on similar in-medium spectral functions) have shown rather good agreement, see e.g. Refs. [1141 130] , 
Ref. [TT5], or Refs. [TrMfTTrJ] . 
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220 MeV (peak) 


250 MeV (peak) 


190 MeV (average) 


T c 


170 MeV 


170 MeV 


175(160/190) MeV 


T fo 


130 MeV 


130 MeV 


120(135) MeV 


spatial 


Glauber (initial) 


Woods-Saxon 


isotropic 




~0.5-0.55 


0.57 


0.53 


v(r) 


approx. linear 


Pt oc \fr 


vt oc r 


tfb 


~ 8-9 fm/c 


~ 7.5 fm/c 


~ 6.5 fm/c 


QGP-EoS 


massless (Nf=3) 


quasi-particle model 


massless (iVy=2.3) 


HG-EoS 


lowest SU(3) multipl. 


m B ,M < 2, 1.5 GeV 


m B , M < 2, 1.7 GeV 


s/qb 


42 


26(?) 


27 


/4° 





? (/0) 


80(35) MeV 


EM rates 


chiral virial 


empirical scatt. ampl. 


hadronic many-body 



Table 1: Fireball parameters employed in the calculations of dilepton spectra in 
In(158 AGeV)-In collisions in Refs. J77?|[II3/ (DZ), {JMUMI (RR) and JM, WUb^ (HR). 

to NA60 data, is compiled in Tab. [T] (see also Ref. |125| ). The overall range of the 
underlying parameters and assumptions is rather similar. This is not a coincidence but a 
consequence of constraints from measured hadron spectra at SPS energies |12H 11221 H23| 
[124] which all of the three models have been subjected to. There are, however, noticeable 
differences. E.g., all approaches operate with a for SPS energies "canonical" formation 
time of Tq=1 fm/c, but the initial peak temperature in Ref. |120j is about 15% larger than 
in Ref. [11 7j (e.g., due to differences in the underlying QGP EoS). Averaging over the initial 
spatial density profile typically leads to a 15% smaller average temperature [126j : thus, 
T ~190 MeV in Ref. [36] is quite consistent with T max ~220 MeV in Ref. pIT]. The slightly 
larger expansion velocity in Ref. |120j (surface velocity t| o =0.57 at thermal freezeout), 
together with its square-root radial profile, imply larger boost factors in the qt spectra 
which becomes significant at high momenta. In this approach preliminary NA60 pion 
spectra in semicentral In-In are saturated by thermal emission over the entire measured 
range up to pr~3 GeV. Alternatively, in Ref. [106] . based on an analysis of pion spectra in 
Pb-Au and S-Au collisions at SPS, the thermal component was found to account for the 
pion yields only up to pt— 1 GeV, requiring the introduction of a "primordial" component 
associated with initial hard scattering of the incoming nucleons. This interpretation is 
supported by the observation that the pion spectra for pr>2 GeV essentially scale with 
the number of binary N-N collisions (iV co ii), indicating that the hard component dominates 
the spectra at these momenta. The preliminary NA60 pion spectra are also well predicted 
in this approach, with a crossing of thermal and hard components at pr— 1.2 GeV. As 
discussed in Sec. II .2|. a valuable indicator of the degree of thermalization is the elliptic flow, 
V2(pt)- At SPS energies, ideal hydrodynamics overpredicts this quantity even at low px by 
about 30-50% (possibly due to neglecting effects of finite viscosity, in connection with initial 
temperatures in the vicinity of T c where the EoS is presumably rather soft). Moreover, the 
experimental V2{pt) i n semicentral Pb-Au levels off at pr=l-5-2 GeV [127^1128"] . indicative 
for a transition to a kinetic regime, while hydrodynamic results keep rising, overpredicting 
^2(pt=2 GeV) by about a factor of ~2. 

Concerning effective chemical potentials for pions (and other stable particles) between 
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chemical and thermal freezeout, their main effect is a faster cooling in the evolution of 
T(t) (the equation of state, P(e), is largely unaffected) [J29l [J30] . E.g., for / u 7r =80MeV 
at T=120MeV [3B], the pion density, g n (T,p n ), is enhanced by a factor of ~e 80 / 120 ~2 
relative to /x 7r =0, and as large as ^ 7r (T=150 MeV,0). Therefore, thermal freezeout at 
(T fo ,^°)=(130,0) MeV [HZ] corresponds to a smaller pion density (and thus larger volume) 
than at (T fo ,^°)=(120,80) MeV [35], consistent with the longer lifetime in Ref. [TT7] . 

Implications of varying critical and chemical-freezeout temperatures for dilepton spec- 
tra have been studied in Ref. |106| . The value of T c affects the relative partition of QGP and 
hadronic emission, especially at masses M>1 GeV where the Boltzmann factor augments 
the sensitivity to earlier phases and the hadronic rates are not enhanced by resonances. 
However, if hadronic and QGP rates are "dual" around T c , this distinction is largely aca- 
demic. Smaller T c h's lead to smaller /Vs in the subsequent hadronic evolution, and thus 
higher kinetic freezeout temperatures, e.g., (Tf Q ,^)=(135,35) MeV for T c h=160 MeV. 

Experimentally measured dilepton spectra contain sources other than thermal radia- 
tion represented by Eq. (|4ip . A systematic evaluation of these sources has recently been 
conducted in Ref. [106] . in terms of (i) final-state decays and (ii) primordial sources. 

Dilepton decays of long-lived hadrons (mostly r], ?/, to and <fi mesons) after thermal 
freezeout, commonly referred to as "hadron decay cocktail", are usually based on chemical 
freezeout for their abundance and thermal freezeout for their pt spectra. The situation is 
more involved for the p-meson, since its continuous regeneration implies relative chemical 
equilibrium with pions until thermal freezeout (to a certain extent this may also apply to 
00 and 4>) . In addition, its short lifetime is not well separated from the typical duration of 
the freezeout process. In Ref. [30], the final generation of p decays has been approximated 
by an extra 1 fm/c of fireball lifetime. However, as has been clarified in Refs. |131[ 1106] . 
when treating the final generation of p's as a cocktail decay, the time dilation of the moving 
p's generates a hardening of its q t spectrum by a factor "f t =Mt/M (Mf=M 2 +q 2 ). The 
resulting spectrum recovers the standard Cooper-Frye [132] description for freezeout at a 
fixed time in the laboratory frame (cf. also Refs. [13311134] ). This, in turn, implies that the 
apparent temperature of the radiation formula (|4ip is smaller than the actual temperature 
figuring into the Boltzmann factor (independent of flow effects) by about ^10%l 10 l 

In analogy to the pion-pr spectra discussed above, the p spectra are expected to have 
a primordial component (emanating from hard N-N collisions) which does not equilibrate 
with the medium. Such a component has been introduced in Refs. [1311 1106] including a 
schematic treatment for Cronin effect and jet-quenching as inferred from pion spectra in 
S-Au and Pb-Au collisions at SPS [135| (also note that, at high q t , this component scales 
with N co \i, rather than A r part as for the (low-pr) cocktail). 

Another primordial dilepton source is the well-known Drell-Yan (DY) process, i.e., 
quark- antiquark annihilation in binary N-N collisions. To leading order qq — ► e + e~ is 
0{oP s a 2 ) and can be reliably calculated in perturbation theory at sufficiently large masses, 
M>2 GeV, utilizing parton distribution functions as input |136j . A finite pair momentum, 
qt>0, can be generated by intrinsic parton kt and at next-to-leading order (NLO) (the 
latter is the dominant effect). The extrapolation of DY to small masses is problematic, 
but at a scale of qt—2 GeV its contribution to dilepton spectra at SPS is potentially sizable. 
In Ref. [106J it has been suggested to estimate the spectrum of slightly virtual DY pairs, 
i.e., for M 2 <C q 2 , by an extrapolation of a finite-^ DY expression to zero mass and 
constrain the resulting photon spectrum by measured photon spectra in p-A collisions. 

10 The time dilation factor for p decays in the thermal radiation formula is compensated by the same 
time dilation in p formation, as a consequence of detailed balance. 
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Semileptonic final-state decays of correlated of D and D mesons (i.e., corresponding to 
an associately produced cc pair) lead to an irreducible dilepton signal. The pertinent mass 
spectrum is, in fact, sensitive to reinteractions of the charm quarks and/or hadrons in the 
medium. At SPS energies, the relevance of this effect for correlated DD decays is currently 
an open question [137] . Theoretical calculations discussed in the following are employing 
iV co ^-extrapolated spectra from p-p collisions (based on PYTHIA [139] simulations) as 
provided by the NA60 collaboration [35]. 

4.2.2 CERN-SPS I: NA60 

In section we discuss several calculations of p + /i _ spectra in semicentral In(158AGeV)- 
In collisions as measured by NA60. The excellent mass resolution and statistics of the 
data allowed for a subtraction of the hadronic cocktail (excluding p and DD decays), 
resulting in the so-called "excess spectra" (in more recent, acceptance-corrected, NA60 
spectra [138] correlated DD decays are also subtracted, with some caveat as to their 
medium modifications, as mentioned above). 

Theoretical predictions [20] of the low-mass excess spectra utilizing the in-medium p 
spectral function of Ref. [30] (cf. Sec. I3,4.2j) showed good agreement with the first data 
of NA60 [35]. More complete calculations including QGP radiation (as in Ref. [140] but 
with hard-thermal loop resummed rates [78]), in-medium u> and (/) decays [105], 4-7r-like 
annihilation (relevant at intermediate mass) [36], as well as primordial p and Drell-Yan 
(DY) contributions (relevant at high qt), are summarized in Fig. 1171 In connection with a 
slight update of the fireball model (larger acceleration implying smaller lifetime), the re- 
sulting description of the NA60 invariant-mass spectra is quite satisfactory over the entire 
range, including projections onto low (gt<0.5 GeV) and high (q t >l.O GeV). In-medium u 
and (j) contributions are rather localized in mass, while QGP and DY radiation are at 
the 10-15% level at masses below 1 GeV. The in-medium plus freezeout (FO) p contri- 
butions [20] remain the dominant source confirming the notion that the NA60 low-mass 
data probe the in-medium p spectral function. This is also borne out of the acceptance- 
corrected <7£-spectra where, for M<lGeV the p contribution prevails up to momenta of 
qt— 1 GeVcj At qt>l GeV DY and primordial p- mesons become an increasingly important 
source, but the data for M=0.4-0.6 GeV and 0.6-0.9 GeV cannot be fully accounted for. 
These discrepancies are less pronounced for central In-In collisions, and may possibly be 
resolved by a stronger transverse expansion within the constraints of the hadronic spectra 
(we return to this question below). At masses M=l-1.4GeV, the most significant sources 
are hadronic emission from multi-pion states (e.g., ir-ai, p-p or ir-u annihilation), QGP 
and DD decays. The hadronic contribution is significantly enhanced (by maximally a 
factor of ~2 around M~l GeV) due to the effects of chiral mixing [36] (recall Sec. I3.2.ip . 
which currently cannot be discriminated by the data. The qt spectra for M=l-1.4 GeV 
are well described over the entire momentum range. 

Fig. [18] summarizes the results of hydrodynamic calculations |117[ 1118] based on 
hadronic emission rates within the chiral-reduction approach |74] [75| [76] (cf. Sec. 13.2 . 1 1) . 
freezeout p mesons (including the proper 7 factor relative to thermal radiation) and per- 
turbative qq annihilation in the QGP (pQGP). The overall structure of the NA60 mass 
spectrum is roughly reproduced (cf. upper left panel of Fig. [18]), but the p resonance figur- 

n The experimental (/(-spectra in Figs. 1171 1181 119l are not absolutely normalized; the theoretical qt-spectra 
in Fig. [17] are normalized using the M-spectra at low qt; however, whereas the experimental M-spectra for 
q t >l GeV are reasonably reproduced, the q t -spectra for M=0.4-0.6 GeV are underestimated for q t >l GeV. 
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Figure 17: Dimuon invariant-mass (left panels) and acceptance- corrected transverse- 
momentum (right panels) spectra in semi-central In(158 AGeV)-In collisions. Calcula- 
tions 1106^ for thermal emission utilizing in-medium p, oj and 4> spectral functions based 
on hadronic-many body theory \3Ch, \10Sj . An annihilation including chiral mixing [36] and 
QGP emission, supplemented by non-thermal sources (Drell-Yan annihilation, primordial 
and freeze-out p-meson, open-charm decays), are compared to NA60 data jl35\, \14A \14^j - 



ing into the EM correlator lacks significant in-medium broadening, despite the reduction 
in peak strength due to the mixing effect (the agreement improves for semiperipheral and 
peripheral collisions |117j ). The freezeout-p contribution compares quite well with the one 
in the upper left panel of Fig. [17] which includes a broadening but also occurs at higher 
pion density (recall the discussion in Sec. 14.2. lj) . The level of the pQGP contribution is 
very similar to the fireball model of Refs. [361 1106] in Fig. [T71 As in Refs. [36] 1106] the 
hadronic contribution at M>lGeV is based on a fit to the EM correlator in vacuum, 
but the mixing effect is less pronounced in the virial scheme, leading to a slightly smaller 
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Figure 18: NA60 dimuon spectra ]35\ 141\\14^ compared to calculations \117Ml3\j employ- 
ing thermal rates from the chiral reduction approach fWj an d perturbative qq annihilation, 
folded over a hydrodynamic expansion for semicentral In(158 AGeV)-In collisions, sup- 
plemented with free p-meson decays after thermal freezeout. Upper left: M -spectra; other 
panels: qt spectra in three mass bins. 



contribution in the dilepton spectrum (possibly also due to the absence of pion chemical 
potentials). The shapes of the qt spectra (local slopes) of all 3 contributions displayed in 
Fig. [18] (in-medium hadronic, freezeout p and QGP) agree well with the fireball calcula- 
tions of Refs. [3611106] as demonstrated in a direct comparison in Ref. [141] . This suggests 
good consistency of the fireball and hydrodynamic evolution. 

In Ref. |120j a thermal fireball expansion (cf. Tab. [H middle column) has been applied 
to compute dimuon spectra utilizing in-medium p and u spectral functions (based on 
empirical scattering amplitudes on pions and nucleons, recall solid lines in Fig. [T2|) [99J, 
vacuum 4-7T annihilation (with both charged |144| and neutral pions), as well as QGP rates 
based on the quasiparticle model of Ref. [143] , cf. Fig. [T9l The overall shape and magnitude 
of the mass spectra is rather well reproduced, except close to the dimuon threshold where 
the importance of baryon effects is apparently underestimated (the underlying p spectral 
function at T=150MeV shows little variation between baryon densities of qb=0-5qo and 
Qo [120| ). The qt spectra can be reasonably well described without contributions from 
DY or primordial p's, which differs from the hydrodynamic (DZ) |117[ I118| and HR- 
fireball [36l EDS] results, cf. Ref. [TH] (recall that the RR fireball model describes NA60 
pion spectra over the entire px range by thermal emission) ; part of this discrepancy is due 
to the slightly larger expansion velocity and the square-root radial profile of the transverse 
rapidity, cf. Tab. [H Another significant difference concerns the magnitude of the QGP 
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Figure 19: NA60 dimuon spectra ]35\ compared to calculations \120$ employing 

thermal rates with in-medium p and u spectral functions Wttf and free 4tt annihilation 
in hadronic matter, as well as qq annihilation in the QGP \143^ , folded over a thermal 
fireball expansion for semicentral In(158 AGeV)-In collisions, supplemented with free p- 
meson decays after thermal freezeout. Left panel: M -spectra for all qt (bottom) and for 
qt>l GeV (top); right panel: qt spectra in three mass bins. 



contribution, which is by a factor of >2 larger in Ref. [120j than in Refs. \36 \ HIT] . Part 
of this discrepancy is due to the quasiparticle QGP EoS employed in Ref. [120] , which 
entails larger temperatures (including To) at given fireball volume. It is also related to the 
prevailing role of QGP radiation for M>1 GeV. 

The sensitivity of the NA60 data to the critical and chemical freezeout temperatures 
has been elaborated in Ref. [106| . by varying T c from 160-190 MeV and T c h from 160- 
175 MeV (keeping the fireball expansion parameters fixed), representing current uncer- 
tainties in lattice QCD [HE] and thermal model fits With "quark-hadron" duality 




Figure 20: Sensitivity of dilepton spectra to critical and chemical- freezeout temper- 
ature llOby . Left panel: QGP and hadronic (4ir) radiation at intermediate mass 
for T c =T ch =175 MeV (EoS-A), T c =T ch =160 MeV (EoS-B) and (T c ,T ch )=(190,160) MeV 
(EoS-C); right panel: slope parameters in qt- spectra for EoS-B and EoS-C (lower curves) 
and for 15% increased fireball expansion (EoS-B + and EoS-C^ , where the latter addi- 
tionally includes u t-channel exchange in irp — ► 7re + e~ reactions; all without Drell-Yan 
contribution ); data are from Refs. \141\ \14%V - 
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Figure 21: Comparison of acceptance- corrected excess dimuon mass spectra (for 
0.2<q t /GeV<04 (left), q t >0.2 GeV (middle) and 1.8<q t /GeV<2.0 (right)) in minimum- 
bias In(158 AGeV)-In collisions \145\ \13l\\138^ to model predictions for semicentral In-In 
of RH (EoS-A ) fTOby . RR \TWj and ZD JiTf . normalized to the average N ch of the data. 



in the thermal dilepton rates [30] in this temperature regime (at all masses, cf. Fig. [T4~]) . 
the invariant-mass spectra turn out to be remarkably insensitive to these variations [106] 
(duality of the QGP and hadronic emission for M < 1.5 GeV close to T c is not realized in 
the rates underlying the calculations of Refs. \117\ 1120"] ). However, the partition of QGP 
and hadronic (4ir) emission at intermediate masses changes appreciably from hadron-gas 
dominated spectra for T c >175 MeV to QGP dominated ones for T c =160 MeV, cf. left panel 
of Fig. [20l In the latter case, the smaller value for 7^=160 MeV implies smaller chemical 
potentials in the hadronic phase. This is part of the reason for the reduction in hadronic 
emission, but also leads to a larger freezeout temperature by about 15 MeV (recall the 
discussion in Sec. 14.2.11 and right column in Tab. [T]). This, in turn, helps in the descrip- 
tion of the transverse-momentum spectra at qt>lGeV. However, an additional increase 
in the transverse fireball acceleration by 15% seems to be required to achieve quantitative 
agreement with the effective slope parameters as displayed in the right panel of Fig. [20] It 
remains to be checked whether this can be consistent with a more complete set of hadronic 
spectra in In(158 AGeV)-In collisions. 

A comparison the three model calculations discussed above to acceptance-corrected 
mass spectra in minimum-bias In(158 AGeV)-In j!45] in Fig. [21] reiterates the importance 
of baryon-driven medium effects [106] at low M and low qt, as well as the lack of high-gj 
yield in the p-mass region and below for Refs. [118] and [106] with EoS-A. The latter 
improves when increasing the fireball expansion as in the right panel of Fig. [20] Also note 
that comparing minimum-bias data to calculations at an average iV c h underestimates the 
theoretical contributions which scale with A^on (DY and primordial p's). 



4.2.3 CERN-SPS II: CERES/NA45 and WA98 

The refinements in the analysis of the NA60 dimuon spectra (fireball evolution and addi- 
tional sources) have been rechecked against existing and updated EM data at the SPS. 

The updated calculations of Ref. [106] agree well with the combined '95/'96 CERES 
dielectron data (left panel of Fig. [22]). For the 2000 data (right panel of Fig. [22]), the 
cocktail-subtracted excess spectra in central Pb-Au corroborate the main findings of the 
NA60 data, i.e., a quantitative agreement with the in-medium p of Ref. [30J and the 
predominance of baryon effects. The longer lifetime of the fireball in central Pb-Au (factor 
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Figure 22: CERES/NA45 data for e + e~ invariant-mass spectra in Pb(158 AGeV)-Au 
collisions at the SPS \32\ \37$ . The left panel corresponds to semicentral collisions including 
the contribution from long-lived hadron decays after freeze- out ("cocktail", dash- dotted 
line), while the right panel shows "excess" spectra for central collisions. The theoretical 
predictions are based on in-medium p spectral functions \3(tf supplemented by u> and <f> 
decays, as well as Drell-Yan and 4^ annihilation 1106\/ (as for the NA60 data, see Fig.\17\). 

~2 relative to In-In) reduces the uncertainties due to p-meson cocktail contributions. In 
addition, dielectrons enable access to very low masses, where the '00 CERES data may 
bear a first hint of a large enhancement as predicted by hadronic many-body theory. 

As emphasized in Sec. 13. 5] (very) low-mass dilepton rates are intimately related to 
thermal photon spectra. In Ref. [108] the in-medium p spectral function of Ref. [30] 
has been carried to the photon point and convoluted over the same fireball expansion 
as before; when supplemented with ^-channel exchange reactions, QGP emission and pri- 
mordial (hard) photons constrained by p-A data, the resulting qt spectra are consistent 
with WA98 photon spectra, see lower right panel of Fig. [23j the updated fireball evolution 
barely affects the total spectra. The contributions from the lightlike p are prevalent up 
to qt—lGeV (cf. Fig. [T5j) . after which t-channel processes takes over. Primordial pho- 
tons outshine the combined thermal yield (hadronic+QGP) for <7i>2 GeV. This is nicely 
consistent with the calculations of Ref. |149j . see lower left panel in Fig. [23l In earlier 
calculations of Refs. [14811126] the thermal yield is significantly larger, due to an increased 
QGP contribution caused by a short formation time of ro=0.2-0.33 fm/c with associated 
peak temperatures of up to To=335MeV (for ro=lfm/c [1081 1149] average initial temper- 
atures are slightly above Tb=200MeV). Even for this upper estimate of QGP emission^!, 
the latter is smaller than the hadronic one for momenta <?t<1.5-2 GeV, and the pQCD 
photons are at the ~40% level of the combined thermal contribution at qt— 2 GeV. 

Similar conclusions arise from theoretical analyses [1401 1150] 1151] of intermediate- 
mass dimuon spectra (1.5>M„„/GeV<3) in Pb-Pb collisions at SPS [15 2 j : unless the 
initial temperature significantly exceeds To=250 MeV, the thermal contribution falls below 
primoridial sources (DY) at masses and transverse momenta beyond M, qt~ 1.5-2 GeV. 
This is fully confirmed by the recent NA60 intermediate-mass dilepton spectra [137] . 

12 At SPS energy, with a Lorentz contraction of 7~9 for the incoming nuclei, the time for full nuclear 
overlap is ca. 0.8 fm/c. 
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Figure 23: Direct photon spectra in central Pb(158 AGeV)-Pb collisions at the SPS as 
measured by the WA98 collaboration 1146[ IM^i - The upper panels contain theoretical 
calculations \14$\ \120j with a moderate primordial pQCD component and a hot initial 
state (tq=0.2-0.33 fm/c; in the upper left panel the lower curves are without primordial 
emission), while the calculations in the lower panels 1149[\103^ contain primordial emission 
with Cronin effect and a larger thermalization time (Vq— 1 fm/c). 



4.2.4 Future Dilepton Measurements 

Dilepton programs will be pursued with high priority over a wide range of collision energies. 
The large enhancement observed in a low-energy (40 AGeV) run at SPS [38] is in line with 
the prediction of hadronic-many body theory that medium effects caused by baryons play 
a leading role [70]. This trend continues down to much lower bombarding energies of 1- 
2 AGeV. However, at these energies recent transport calculations suggest that the low-mass 
enhancement, which could not be explained by hadronic in- medium effects [27], is related 
to primordial iV-iV Bremsstrahlung |153j . as well as A — > Ne + e~ and n Dalitz decays [1541 
155J. A better sensitivity to medium effects appears to be in the p-uj mass region, where 
the (lack of) yield indicates a strong broadening of the vector resonances [1561 1154] 1157] . 

At the high-energy frontier, first RHIC data [42] find a large e + e~ signal especially in 
the mass region around M~0.3 GeV. The excess is concentrated at low qt and in central 
collisions, and cannot be explained by current in-medium spectral functions. It is tempting 
to speculate that the excess is caused by the formation of a disoriented chiral condensate 
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(DCC), as pion-DCC annihilation shares the above features |158| , However, the magnitude 
of this dilepton source cannot easily compete with hadronic medium effects, unless the 
DCC domains are rather large and abundant. In this case, footprints of the DCC should 
be visible in other observables (e.g., 7r° and px fluctuations). Precision measurements 
within the RHIC-II program will be of crucial importance here |43j . 

Finally, dilepton data will play a critical role in the CBM experiment at the future 
GSI facility (FAIR). In the planned energy regime, £^=10-40 AGeV, one envisages the 
largest nuclear compression and thus maximal baryon density, ideally suited to scrutinize 
the current understanding of medium effects. An extra benefit could be the occurence of 
a critical point or a true mixed phase at a first-order transition, with extended fireball 
lifetimes further enhancing the dilepton signal. 

4.3 Critical Appraisal 

In this section we evaluate the current status of determining the in-medium vector spectral 
functions (focusing on the p meson) and the implications for chiral restoration. 

Calculations of p-meson spectral functions based on effective chiral Lagrangians cou- 
pled with many-body techniques agree on a strong broadening with small (positive) mass 
shifts. At normal nuclear matter density, one finds an increase in width of Ar p ~250 MeV 
with an estimated error of ~30%, i.e., the vacuum width almost triples. The question 
whether the parameters in the effective Lagrangian are subject to in-medium changes 
requires further input. In the vector-manifestation scenario, reduced bare masses and 
coupling constants are inferred from a matching of the correlators to an operator prod- 
uct expansion (OPE) at spacelike momenta governed by the in-medium reduction of the 
condensates. However, it turns out that, within current uncertainties, the softening of the 
p spectral function as imposed by the OPE at nuclear matter density is fully accounted 
for by the broadening due to hadronic many-body effects. More accurate tests of this 
assertion, especially at higher densities/temperatures, will require a more precise determi- 
nation of the in-medium condensates on the OPE side of the QCD sum rule. The predicted 
broadening is supported by several recent experiments where dilepton spectra have been 
measured with impressive precision: at JLAB, photoproduction data off mid-size nuclei 
find a p broadening of Ar^ uc ~70-100 MeV without significant mass shift, consistent with 
many-body effects at about half nuclear density and 3-momenta of ~l-2 GeV. The NA60 
dilepton spectra in central In-In collisions exhibit an average p width of r™ c ~400MeV 
i.e., an additional broadening of Af^ IC ~250 MeV. Typical kinetic freeze-out conditions 
at SPS energies are (p§,T fo )~(0.3yo ,120MeV). With initial temperatures of T ~200MeV 
(as suggested by "effective" slope parameters in the qx spectra for M>lGeV, as well as 
direct photon spectra in Pb-Pb), the average p width thus reflects the medium at an aver- 
age temperature of T~ 150 MeV (the growing fireball 3- volume "biases" low-mass dilepton 
radiation to more dilute stages). This implies that the p width approaches its mass when 
the system moves toward the (pseudo) critical temperature, T p (T— >T C )— > m p , i.e., the 
resonance "melts" (see also Ref. |159j ). Inspection of the theoretical predictions for the 
width of the p as extracted from the vector spectral function corroborates this conclusion, 
cf. Fig. [24l The circumstantial "duality" of hadronic and partonic EM emission rates 
close to T c lends robustness to the pertinent predictions for dilepton spectra in heavy-ion 
collisions as they become independent on details of the evolution model, in particular of 
the treatment of the phase transition region. The excess radiation at intermediate mass, 
with its rather soft emission characteristics in q t , as well as direct photon spectra, further 
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Figure 24: Temperature dependence of the p width as evaluated in hadronic many-body 
theory 130^ at fixed entropy per baryon representative for heavy-ion collisions at full SPS 
energy (dashed and solid line: with and without meson chemical potentials). For the dash- 
dotted line, medium effects due to baryons are switched off. 

consolidate the origin of thermal radiation from around T c . 

A no less challenging task is to connect the above findings to the chiral transition. It has 
recently been argued that the statistical operator consisting of a (free) hadron resonance 
gas (HRG) is capable of reproducing several features of lattice QCD computations of the 
equation of state (EoS) until rather close to the (pseudo-) critical temperature [1601 1161] . 
Beyond T c the lattice EoS levels off, reflecting quark-gluon degrees of freedom, while the 
HRG EoS diverges (Hagedorn catastrophe). The p melting offers a microscopic explana- 
tion for this transition: under moderate conditions, the interacting HRG physics drives 
the p broadening to an extent which justifies the use of well-defined quasi-particle states 
in the statistical operator. With further increasing temperature and density, resonance 
overlap in the p spectral function drives it to a continuum shape with a strength resem- 
bling a weakly interacting qq pair, i.e., the resonance strength in the statistical operator 
converts into partonic strength. The phenomenon of overlapping resonances merging into 
a perturbative qq continuum is, of course, well known from the e + e~ annihilation cross 
section into hadrons above M~1.5 GeV. It is suggestive that the thermal medium provides 
the necessary phase space for low-mass resonances which, via their mutual "mixing" in 
different hadronic correlators, "restore" quark-hadron duality down to M^0, implying 
chiral restoration. To quantify this picture the evaluation of chiral order parameters is 
mandatory. It is tempting to speculate that the rather sharp increase of the p width 
close to the expected critical temperature (especially in chemical equilibrium as realized 
in lattice QCD, represented by the solid line in Fig. [24"j) is signaling the chiral transition. 
QCD sum rules remain a valuable tool if the T (and pb) dependence of the quark and 
gluon condensates can be made more precise. Ideally, the latter are determined from first- 
principle lattice QCD calculations. Possibly the most promising approach, which has been 
little exploited thus far, are chiral (or Weinberg) sum rules. Their use hinges on the in- 
medium axialvector spectral function. The latter is much more difficult to constrain due to 
a principal lack of experimental information, encoded in either 3-pion or n-'y final states. 
This stipulates the importance of calculating the axialvector correlator in chiral models. 
In connection with a realistic vector correlator and lattice-QCD input on the in-medium 
condensates, the explicit realization of chiral restoration can be investigated. First efforts 
in this direction have been undertaken |162^ I163| 1164^ [T65] . but a full treatment including 
quantitative V and A spectral functions, even in the vacuum, is currently lacking. 
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5 Conclusions 

Medium modifications of hadronic spectral functions play a key role in the diagnosis of 
hot/dense strongly interacting matter and its condensate structure. Experimentally, the 
most promising approach is dilepton spectroscopy which directly probes the vector spectral 
function of the hadronic medium. For the p-meson, which dominates the low-mass vector 
channel, effective hadronic theories largely agree on a strong broadening of the resonance, 
with little mass shift. Baryon effects prevail over those induced by mesons, and the 
predicted modifications in cold nuclear matter are compatible with QCD sum rules at 
finite density Intense experimental efforts over the last ~15 years have culminated to 
a new level of precision which broadly confirms the theoretical expectations: production 
experiments off ground-state nuclei find an increase of the p width by ~80 MeV, while the 
effect in heavy-ion collisions at the SPS is by a factor of ~3 larger. Part of this difference 
is due to the access to the low-momentum regime in the heavy-ion measurements. It 
is therefore highly desirable to push the sensitivity of the nuclear experiments to low 3- 
momenta where significantly larger medium effects are predicted. The average p width 
extracted in heavy-ion collisions suggests that the p resonance "melts" close to the expected 
phase boundary, in agreement with extrapolations of hadronic models. This is a first 
explicit evidence that melting resonances are involved in the transition from hadronic to 
quark degrees of freedom. Modern quark-model calculations could provide complementary 
insights when approaching T c from above. Unquenched lattice QCD computations of the 
vector correlator would undoubtedly set valuable benchmarks and possibly shed light 
on the conjecture that the width is connected to order parameters of chiral symmetry 
restoration. In addition, information on quark condensates and pion decay constant (s) 
below T c can be connected to hadronic vector and axialvector spectral functions utilizing 
Weinberg sum rules. The synergy of hadronic and quark models with first-principle lattice 
QCD computations, augmented by quantitative applications to experiment at current and 
future facilities, opens exciting perspectives to improve our knowledge about the chiral 
transition in hot /dense QCD matter and the generation of luminous mass in the Universe. 
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